1/1 


RcMUM  719  THE  EFFECTS  OF  OPTIC*.  SMOOTHING  TECHNIQUES  ON 

FILAHENTRTION  IN  LASER  PLASHASOJ)  NAVAL  RESEARCH  LAt 
HASHINOTON  DC  A  J  SCHHITT  29  DEC  97  NRL-HR-C142 

r/a  29/9 


UNCLASSIFIED 


Naval  Research  Laboratory 


Washington,  DC  20375-5000 


OTIC  File  copy 


NRL  Memorandum  Report  6142 


The  Effects  of  Optical  Smoothing  Techniques  on 
Filamentation  in  Laser  Plasmas 


Andrew  J.  Schmitt 


Laser  Plasma  Branch 
Plasma  Phxsics  Division 


December  29,  1987 


DTIO 

PIHLECTEI 


Approved  for  public  release,  distribution  unlimited 


- .  <.  r.  /. 


SECURITY  CLASSIFICATION  OF  *HiS  RAGE 


la  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2a  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

NRL  Memorandum  Report  6142 


6a.  NAME  OF  PERFORMING  ORGANIZATION 


REPORT  DOCUMENTATION  PAGE 


lb  RESTRICTIVE  MARKINGS 


Form  Approved 
OMB  No  0704-0188 


3  DISTRIBUTION/AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 
unlimited . 


5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


Naval  Research  Laboratory 


6c  ADDRESS  (C/ty,  State,  and  ZIP  Code) 
Washington,  DC  20375-5000 


6b  OFFICE  SYMBOL  I  7a  NAME  OF  MONITORING  ORGANIZATION 
(If  applicable)  | 


Code  4730 


7b  ADDRESS  (C/ty,  State,  and  ZIP  Code) 


8a.  NAME  OF  FUNOING/ SPONSORING 
ORGANIZATION 

Department  of  Energy 


8c  ADDRESS  (City,  State,  and  ZIP  Code) 

Washington,  DC  20545 


8b.  OFFICE  SYMBOL  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

(If  applicable) 


to  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 


WORK  UNIT 
ACCESSION  NO 


11.  TITLE  (Include  Security  Clarification) 

The  Effects  of  Optical  Smoothing  Techniques  on  Filamentation  in  Laser  Plasmas 


12.  PERSONAL  AUTHOR'S) 
Schmitt,  Andrew  J. 


13a.  TYPE  OF  REPORT 
Interim 


16  SUPPLEMENTARY  NOTATION 


13b  TIME  COVERED  14  DATE  OF  REPORT  [Year.  Month  Day)  IS  PAGE  COuU’ 

from  10/86  to  10/87  1987  December  29  oj 


COSATI  CODES 


GROUP  SUB-GROUP 


*8  SUBJECT  TERMS  {Continue  on  reverse  it  necessary  j  no  identity  Dy  o/o  ck  n>jmoeri 

Laser  plasma  interaction 
Inertial  confinement  fusion 


'9  ABSTRACT  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

The  effect  of  the  induced-spatial-incoherence  (IS!)  and  the  random-phase-screen  (RPS)  optical  smoothing  techniques 
on  the  filamentation  instability  in  laser  plasmas  has  been  investigated  numerically  and  analytically.  A  two-dimensional 
time-dependent  laser-plasma  propagation  code,  including  both  ponderomotive  and  thermal-conduction  dominated  tiiamenta- 
tion,  is  used  to  simulate  the  laser-plasma  interaction.  The  results  of  these  simulations  are  compared  to  the  predictions  of  a 
simple  theory  that  describes  the  filamentation  of  both  coherent  light  and  spatially  and  temporally  incoherent  light.  It  is 
shown  that  filaments  driven  by  the  thermal  mechanisms  tend  to  cluster  together  and  produce  greater  large-scale  nonunifor¬ 
mities  in  the  laser  illumination  than  the  ponderomotively-dnven  filaments.  The  RPS  optical  smoothing  technique  is  tound 
to  reduce  filamentation  only  if  fast  focusing  optics  (F  <  5)  are  used.  The  ISI  smoothing  method  suppresses  filamentation 
for  fast  or  slow  focusing  optics,  and  requires  only  moderate  laser  bandwidth  iAaj/w  =  0. 1  f»).  In  general,  the  iSl  smooth¬ 
ing  method  provides  the  best  suppression  of  filamentation.  Under  common  laboratory  conditions,  filamentation  is  pro¬ 
nounced  at  longer  laser  wavelengths  ( 1 .06  nm  -  0.53  pun),  suggesting  that  current  experiments  may  be  dominated  by  fila¬ 
mentation  effects.  Optical  smoothing  methods  can  suppress,  but  not  eradicate,  filamentation  effects  at  these  wavelengths. 


J0  DlSTRiauT  ON  Ava.lABil.'y  OF 

ab$~Rac  " 

IS  /NCLASSiF'EQ  UNLIMITED  C 

;ame 

□  rc  ae^s 

*.2 a  'J  AME  jt  jr'/'jVBLl 

gal 

Andrew  J.  Schmitt 

DO  Form  1473,  JUN  36 


Ads  PAL  - C ■_  ,RIT  -  . 

UNCLASSIFIED 


’  \r  .  n< -utJP  t  JOP)  I  .  . 

(l')2)  ''17-368 1 


•Vpvtous  jit  ions  j  re  ■  osoiete 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


'  19.  ABSTRACT  (Continued) 

The  optical  smoothing  methods  are  most  effective  in  short*wavelength  (0.25  /im)  laser-driven  plasmas.  At  0.25  tun  laser 
wavelength,  ISI  is  found  to  completely  eliminate  filamentation  effects  in  both  time-averaged  and  instantaneous  intensity  dis¬ 
tribution. 


Accession  For 

~ NTIS  GRA&I 

* 

DTIC  TAB 

rl 

Unannounced 

□ 

,  Justification — 

By - - - 

Distribution/ _ 

Availability  Codes^ 
j  jAvail  and/or 

1 D 1 s  t  I  Special 


m 


QtMl.lTY 

'  :,f  CT-  . 


S3 


I  4**  -»*■*  * .1 


■  - 1  A‘  m  | 


,  I»4  ».l  ».l  >4|V 


:.  *  k.m  Li  Lj 


»■!  «>i  ->4  aI>  .U  .11  . 


£M 

)*rfl 

$ 

Ijl 


CONTENTS 


INTRODUCTION 


GOVERNING  EQUATIONS  AND  NUMERICAL  DESCRIPTION  .  3 

A.  Light  Propagation  .  3 

B.  Plasma  Response  .  4 

ANALYTIC  FORMULATION  OF  FILAMENTATION  .  8 

A.  Basic  theory  of  filamentation  .  8 

B.  ISI  and  filamentation .  9 

ANALYSIS  AND  DISCUSSION  .  11 

A.  Filamentation  in  Homogeneous,  Nonabsorbing  Plasmas: 

Gaussian  Filaments  .  11 

B.  Filamentation  in  Homogeneous,  Nonabsorbing  Plasmas: 

Complex  Laser  Beams  .  13 

i.  Definition  and  Construction  of  Complex  Laser  Beams  .  13 

ii.  Filamentation  of  Generic  Laser  Beams  .  15 

iii.  Filamentation  of  RPS  Laser  Beams  .  16 

iv.  Filamentation  of  ISI  Laser  Beams  .  16 

C.  Filamentation  of  Complex  Laser  Beams  in  Laboratory  Plasmas  .  18 


CONCLUSIONS 


ACKNOWLEDGEMENTS 


APPENDIX  .  23 

A.  1  Ponderomotive  Filamentation  .  23 

A. 2  Thermal  Filamentation  .  24 


REFERENCES 


TABLE  I  .  29 

TABLE  II  .  30 


FIGURES 


‘'SBI&SWKWSfli 


SI 


THE  EFFECTS  OF  OPTICAL  SMOOTHING  TECHNIQUES  ON 
FILAMENT ATION  IN  LASER  PLASMAS 


I.  Introduction 

The  production  of  uniform  laser  illumination  has  been  a  fundamental  concern  in  the  laser-fusion 
community  for  many  years.  The  quality  of  the  laser-plasma  coupling  and  symmetry  of  the  pellet 
implosion  are  dependent  on  the  illumination  uniformity.  Illumination  uniformity  analysis  can  be 
divided  naturally  into  large  and  small  scale  categories.  Large  scale  illumination  uniformity  is 
primarily  determined  by  the  gross  intensity  profile  of  the  incident  laser  beams,  the  inter-beam  energy 
balance,  and  the  geometry  of  the  targeting  on  the  pellet.  This  aspect  of  the  uniformity  problem  has 
been  previously  addressed,  and  the  results  indicate  that  the  gross  uniformity  can  probably  be 
controlled  to  within  tolerable  limits^.  The  detailed  structure  of  the  individual  laser  beam  is 
responsible  for  small  scale  noauniformities,  and  is  more  difficult  to  suppress  or  control.  The  laser- 
plasma  filamentation  instability  is  seeded  by  these  small  scale  nonuniformities  and  tends  to  increase 
the  nonuniformity,  exacerbating  the  problem.  ^ 

The  filamentation  instability  is  a  nonlinear  optical  effect  that  has  been  studied  for  many  years.1' 
Filamentation  is  caused  by  perturbations  or  nonuniformities  in  light  that  produce  local  changes  in  the 
dielectric  constant,  or  index  of  refraction,  of  a  medium.  If  the  change  in  the  dielectric  constant  is 
positive  in  regions  of  higher  intensity,  a  focusing  lens  is  produced  there.  This  increases  the 
perturbation  and  starts  the  instability.  The  instability  will  saturate  when  the  focusing  tendency  of 
the  intensity  hot  spot  is  balanced  by  diffraction,  but  at  this  point  the  filament  intensity  may  be 
orders  of  magnitude  higher  than  the  initial  perturbation. 

In  laser  plasmas, 
dielectric  constant  and 
heating  or  thermal  effec 

the  dielectric  constant  by  expelling  the  plasma  density  from  the  high  intensity  region;  these 
mechanisms  are  active  in  laser-fusion  plasmas.  The  relativistic  mechanism  affects  the  dielectric 
constant  via  the  increase  in  electron  mass  from  the  relativistic  quiver  velocity  in  the  electromagnetic 
wave;  this  is  generally  unimportant  at  the  intensity  ranges  used  in  laser  fusion  applications. 

Intensity  hot  spots  caused  by  filamentation  of  the  laser  can  seriously  degrade  the  laser-plasma 
coupling,  and  may  affect  the  ablation  pressure  uniformity.  The  high  intensity  filaments  can  induce 
other  harmful  laser-plasma  instabilities,  such  as  stimulated  Raman  scattering  or  two  plasmon  decay. 
Since  filamentation  has  one  of  the  lowest  instability  thresholds,  and  produces  conditions  favorable  for 
other  instabilities,  it  effectively  lowers  the  thresholds  and  increases  the  growth  rates  of  these  other 
instabilities.  These  instabilities  degrade  the  coupling  quality  by  producing  superthermal  electrons 
which  can  preheat  the  fuel  and  spoil  the  gain.  The  enhancement  of  secondary  instabilities  also 
obscures  our  understanding  of  the  underlying  physics  of  the  laser-plasma  interaction:  since 
filamentation  is  itself  hard  to  measure  directly,  it  is  even  more  difficult  to  diagnose  the  effect  of 
filamentation  on  other  phenomena.  There  is  also  the  concern  that  the  laser  filaments  may  produce 
ablation  pressure  nonuniformities  which  seed  or  drive  the  Rayleigh-Taylor  instability  in  the  imploding 
pellet. 

To  avoid  the  filamentation  instability,  as  well  as  to  provide  some  control  over  the  gross  laser 
beam  profile,  new  optical  smoothing  techniques  have  been  developed.  One  of  these  is  the  induced 
spatial  incoherence  ( IS  I )  method  ;  another  is  the  random  phase  screen  (RPS)  method. 
Experimental  results  using  these  techniques  are  incomplete  and  are  still  under  investigation.  The  IS  I 
results  to  date  are  favorable  ;  diagnostics  of  most  laser-plasma  instabilities  (e.g.,  2W^,  3W  3.  Raman 
scatter,  SBS  scatter)  show  significant  reductions  when  compared  to  resuits  with  an  unsmoothed  Seam. 

The  RPS  method  has  also  shown  some  reduction  in  3W  / 2  emission  when  used  in  a  spherical 
13  ° 

illumination  geometry. 


there  are  a  variety  of  mechanisms  that  give  rise  to  an  intensity-dependent 
produce  filamentation.  Among  these  are  ponderomotive  force  effects  ’  ,  plasma 
its®,  and  relativistic  effects.^  The  ponderomotive  and  thermal  effects  change 
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The  implementation  of  each  of  these  optical  smoothing  methods  involves  trade-offs.  The  RPS 
technique  requires  fast  focusing  optics  to  work  effectively,  but  can  be  implemented  with  only  minor 
modifications  to  existing  lasers.  ISI  promises  to  be  a  more  robust  smoothing  technique,  but  requires 
a  broadband  laser  (A Ll/U  ~  0.1%).  Glass  lasers  can  be  converted  to  run  broadband  at  1.06 flm  and 
0.53}lm  laser  wavelengths;  for  short  wavelengths,  the  KrF  laser  is  an  excellent  candidate  for  ISI.  ^ 

Evaluation  of  these  optical  smoothing  techniques  requires  information  on  how  the  filamentation 
instability  is  modified  by  the  incident  laser  beam  structure.  Experimental  data  on  filamentation  has 
been  sparse  or  nonexistent  because  of  difficulties  in  controlling  and  diagnosing  both  the  laser  intensity 
structure  and  the  plasma  conditions,  and  theoretical  techniques  are  insufficient  to  handle  the 
complexity  and  nonlinearity  inherent  in  the  problem.  Computational  techniques  are  also  severely 
tested,  especially  when  modelling  ISI  or  RPS.  Both  ISI  and  RPS  require  resolution  of  a  large  range 
of  scalelengths  (typically  from  less  than  one  to  a  few  hundred  laser  wavelengths).  This  requires  large 
computational  meshes,  and  the  resolution  constraints  invalidate  the  use  of  well  developed  ray-tracing 
computational  methods,  which  ignore  diffractive  effects  that  are  important  at  the  small  scales. 

We  will  analyse  the  limitations  of  optical  control  techniques  and  compare  these  smoothing 
methods  with  traditional  (unsmoothed)  high-power  lasers.  The  paper  is  organised  as  follows:  in 
section  II,  we  describe  the  construction  of  a  two  dimensional,  time  dependent,  laser-plasma 
propagation  code  that  includes  both  ponderomotive  and  thermal  filamentation  mechanisms.  A  steady 
state  version  of  this  code  is  also  developed  to  treat  problems  without  inherent  time  dependence. 
Next,  a  general  analytic  perturbation  formalism  of  filamentation  is  presented  in  section  III,  and  is 
extended  to  account  for  filamentation  of  incoherent  ISI  laser  light.  Finally,  section  IV  presents  the 
results  of  the  numerical  codes  for  a  variety  of  laser  beam  profiles,  and  compares  these  results  to  the 
analytic  predictions. 

We  will  show  that  there  is  a  qualitative  difference  between  thermal  and  ponderomotive 
filamentation:  filaments  created  by  the  thermal  mechanism  tend  to  bunch  together  and  cause  greater 
nonuniformity  than  ponderomotive  filaments,  which  interact  less.  This  clustering  can  result  in  a 
different  saturation  mechanism  for  thermal  filamentation.  We  also  find  that  the  RPS  optical 
smoothing  method  is  dependent  upon  the  use  of  fast  optics  (F/#  £5),  and  that  the  ISI  method  is 
capable  of  suppressing  filamentation  effects  and  providing  smooth  time-averaged  intensity  distributions 
using  moderate  laser  bandwidth.  In  spite  of  the  smoothing  effects,  ISI  can  also  produce  noticeably 
enhanced  intensities  in  the  plasma.  In  laboratory  conditions,  however,  ISI  is  shown  to  suppress 
filamentation  effects  more  completely  than  RPS  or  generic  laser  beams.  At  longer  laser  wavelengths 
filamentation  occurs  readily,  and  optical  smoothing  techniques  suppress,  but  do  not  eliminate,  the 
filamentation  tendency.  The  best  results  for  all  optical  methods  are  found  with  short-wavelength 
(\Q~0.25flm)  laser- irradiated  plasmas,  where  high  absorption  helps  to  reduce  filamentation.  In  these 
plasmas,  the  ISI  technique  can  completely  eliminate  filamentation. 
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II.  Governing  Equations  and  Numerical  Description. 

In  this  section,  we  will  develop  the  basis  of  the  numerical  scheme  and  discuss  the  relevant 
physics  involved.  We  first  derive  the  equations  describing  the  laser  light  propagating  in  a  plasma  in 
which  the  dielectric  constant  has  been  perturbed;  then  we  derive  the  equations  describing  the 
perturbation  of  the  dielectric  constant  caused  by  changes  in  the  plasma  density  responding  to  the 
laser  light.  We  will  treat  two  mechanisms  responsible  for  plasma  density  changes.  The  first  is  the 
ponderomotive  force,  which  directly  expels  the  plasma  from  regions  of  high  laser  intensity.  The 
second  is  thermal  conduction  dominated  plasma  heating,  which  creates  temperature  (and  thus  pressure) 
gradients  and  also  forces  the  plasma  ant  of  high  intensity  regions. 

The  numerical  algorithms  that  handle  the  governing  equations  will  also  be  outlined  in  this 
section.  Two  versions  of  this  code  have  been  created.  The  first  is  time-dependent  and  is  used 
primarily  for  the  ISI  calculations.  The  other  model  solves  the  equations  in  a  quasi-steady  state 
approximation,  and  is  used  to  calculate  filamentation  effects  when  the  incident  light  is  stationary  in 
time. 

The  analysis  in  this  paper  uses  a  two-dimensional  (2D)  cartesian  description  of  the  interaction. 
The  spatial  growth  rates  and  amplification  wavenumber-spectrum  are  very  similar  in  two  and  three 
dimensions,  as  we  will  show  in  section  III.  The  2D  numerical  analysis  also  reduces  computational 
memory  and  time  constraints  to  manageable  levels,  and  allows  a  large  region  of  parameter  space  to 
be  covered  by  using  many  simulations.  Cartesian  rather  than  cylindrical  geometry  is  needed  to  study 
filamentation  in  2D,  since  cylindrical  geometry  artificially  favors  on-axis  focusjng  and  cannot  equitably 
treat  the  random- phasing  required  in  optical  smoothing  simulations  ,i<.  The  main  difference 
between  the  2D  calculations  and  three-dimensional  (3D)  calculations  occurs  in  the  peak  intensity 
values:  3D  filaments  can  typically  reach  much  higher  peak  intensities.  However,  rapidly  varying 
intensity  profiles  may  restrain  large  peak  intensities  in  3D:  preliminary  results  of  3D  simulations  of 
ISI  show  that  the  light  energy  distribution  is  comparable  to  2D  simulations.^^ 

A.  Light  Propagation 

The  laser  propagation  and  filamentation  will  be  described  by  the  parabolic  wave  approximation 
to  the  Maxwell  wave  equation.  Starting  with  the  Maxwell  equations,  we  assume  that  the  divergence 
of  the  electric  field  is  small  (ka*Ve/e<<k^)1  and  set  V^E^O.  The  electric  field  is  separated  into 
fast  and  slow  space-time  scales  with  the  substitution:  E(x,z,t)  =  l/2l'J(x,z,t)exp{-iJkodz  +  iw^tj-f-c.c.i. 
$(x,*,t)^is  the  wave  envelope  of  the  electric  field,  kQ(z)  is  the  (real)  laser  wavenumber 
(^o(*)=u0eor'*^'/c  Cor  “  real  part  unperturbed  plasma  dielectric  constant  (6^  =  6  -r  i£^., 

and  C  (z)  =  l-n  (z)/n^,  where  n<-(z)  is  the  unperturbed  plasma  electron  density,  and  n_  is  the  plasma 
critical  density).  (The  subscript  ’o’  appended  to  a  variable  means  that  it  is  evaluated  at  the 
background  or  unperturbed  state.)  If  the  wave  envelope  (tc, 3 , t )  varies  slowly  with  respect  to  the 
laser  wavelength  and  frequency,  the  Maxwell  wave  equation  reduces  to  the  parabolic  wave  equation: ^ 

2 

(2iko  h  -  ’!>  *  -  (-!  t)-<orw>  -  4ko)  *  v-v> 

c 

*>  2  **  * 

where  ej*,x,t)=  1-W“(a, x,t)/Wo-iU*(z,x,t)l/?i(i,x,t)/Uo  is  the  fully  perturbed  plasma  dielectric  constant; 
U^=4Te”a^/ is  the  electron  plasma  frequency,  and  is  the  electron-ion  collision  frequency.  (The 
time  derivative  can  be  formally  eliminated  by  transforming  to  the  frame  moving  with  the 

pulse  group  velocity*  ’*  ,  with  the  variable  substitution  t’  =  t-i/v  .  In  practice,  one  can  simply  ignore 
the  time  derivative  term  if  the  transit  time  of  the  propagation  region  is  much  smaller  than  any 
characteristic  time  for  changes  in  the  dielectric  constant.)  We  normalize  all  spatial  coordinates  by 
the  laser  vacuum  wavelength,  X  ,  and  define  the  transformation  (using  the  spatially  normalized 
coordinates): 

1  Cz  ~ _  I 

V  =  i  exsIpH  dz*  k  (2re  .(2*)  -  3k/as*n  (2.2', 

•  '2J  0  01'  1  > , 

*  q  ■* 

where  <  sx  z ;  J  anl  £  —-J~  lz;V  iz\u"  ;s  the  im.aginarv  part  •:  th<*  uinerturi-ed  si. ism  a 
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inhomogeneity.  Redefining  the  variable  in  the  i  direction  as  f7=/Jdz’/k(z’),  we  arrive  at  the 
parabolic  wave  equation  in  canonical  form: 


(4Ti3n  'V  *  =  4r^e(7,x,t)  i 


(2.3) 


where  S€=€(T],x, t)-€(z)  is  the  change  in  the  plasma  dielectric  constant  induced  by  the  laser  EM  field. 
The  effect  of  the  plasma  on  the  propagating  field  is  concentrated  in  this  term.  For  the  cases  of 
interest  here,  this  induced  change  is  entirely  due  to  changes  in  the  plasma  density  (<?€  = 
(n(s)  n(7,x,t))/nc),  and  is  found  by  solving  the  equations  governing  the  plasma  response  to  the  laser. 

The  parabolic  wave  equation  (2.3)  is  solved  numerically  on  a  mesh  in  x-z (T])  space,  with  time 
treated  as  an  independent  parameter.  The  numerical  algorithm  consists  of  a  split-step  fast-Fourier- 
transform  (FFT)  technique  .  To  accommodate  the  Fourier  transform  technique,  the  computational 
mesh  is  periodic  and  equally  spaced  in  the  transverse  (x)  direction. 

The  periodicity  of  the  mesh  affects  the  manner  in  which  the  initial  laser  field  is  numerically 
constructed.  The  focusing  optics  separate  the  incident  laser  beam  into  many  individual  beamlets, 
which  are  focused  onto  a  target  where  they  overlap  one  another.  In  the  code,  these  overlapped 
beamlets  are  approximated  as  plane  waves  incident  on  the  plasma  at  different  angles;  this  is  a  good 
approximation  in  the  center  of  the  focal  spot.  For  a  high-power  laser  without  optical  smoothing, 
these  beamlets  correspond  to  light  from  different  sections  of  the  beam,  and  they  will  be  be  slightly 
incoherent  with  respect  to  one  another  (depending  upon  the  degree  of  aberration  in  the  beam).  In 
the  case  of  the  1ST  and  RPS  smoothing  techniques,  these  beamlets  are  incoherent  with  respect  to  one 
another,  and  correspond  to  the  light  coming  from  different  echelon  steps  or  phase  shift  regions  of  the 
random  phase  screen. 

The  incident  wave  field  is  constructed  on  the  Fourier  transformed  plasma  mesh  which  has  the 
spacing  Ak^=2T/X  ajt,  where  is  the  length  of  the  mesh  in  the  x  direction.  Each  beamlet  has 

a  unique  wavenumber  k  ;  the  wavenumber  of  the  i-th  plane-wave  beamlet  incident  at  angle  9.^  has  a 
k  component  kx.=kosinfl-.  The  angles  0.  are  given  by  the  relation  0.-tan  1  (d[i — ( /), 
where  d  is  the  width  of  t\ie  echelon  step,  ^hase  shift  area,  or  spatial  coherence  distance;  /  is  the 
focal  length  of  the  lens,  and  i  varies  from  1  to  (N^  is  the  number  of  beamlets  used).  For 
moderate-to-large  F/#  optics,  is  small  and  tanflj-sinfl.-fl-.  The  modes  corresponding  to  the 

individual  beamlets  are  then  assigned  to  the  nearest  point  on  the  k-space  transform  mesh;  the  mesh 
spacing  Ak^^T/X^^  tbe  Fourier  transformed  mesh  corresponds  to  an  spatial  angular  resolution 
Ad?=X^/X  .  The  real  and  imaginary  parts  of  the  electric  field  of  each  beamlet  at  wavenumber  k^. 
are  then  adcled  to  the  value  at  the  assigned  mesh  point.  In  general,  the  electric  field  at  each  k 
mesh  point  at  i=0  is  independently  assigned  a  random  amplitude  or  phase  (or  both)  depending  upon 
the  type  of  beam  being  simulated.  For  ISI  simulations,  the  phases  of  each  beamlet  are  chosen  to  be 
uniformly  random  and  their  amplitudes  are  selected  with  a  Gaussian  probability  distribution.  For 
RPS  and  other  lasers,  we  typically  use  constant  amplitude  but  uniformly  random-phased  modes. 
(More  detail  on  the  construction  of  laser  profiles  is  given  in  section  IV).  The  laser  electric  fieiu 
E(x,.=0,t  )  is  then  found  by  inverse  Fourier  transforming  this  distribution.  For  the  ISI  simulation, 
this  process  is  repeated  again  after  every  interval  t£  to  produce  a  new  randomly  generated  E  field. 

B.  Plasma  Response  t 

The  nonlinear  change  in  the  dielectric  constant,  5e(x,z,t,f  ]f),  is  found  by  solving  for  the 
plasma  density  using  a  one-fluid  plasma  model.  The  calculation  is  greatly  simplified  by  ignoring 
fluid  coupling  to  itself  along  the  direction  of  propagation  of  the  laser  (z  axis).  This  approximation  is 
valid  when  the  plasma  gradients  along  the  z-axis  are  much  smaller  than  gradients  perpendicular  to 
the  i-axia,  and  is  consistent  with  the  slowly-varying-enveiope  approximation  used  previously  for  the 
laser  electric  field.  The  approximation  also  ignores  fluid  flow  along  the  laser  axis  (which  is  typically 
supersonic  in  the  underdense  regions  of  laser-fusion  plasmas”'*).  The  effect  if  How  on  filamentation 
has  been  considered  elsewhere"^  and  can  be  ignored  if  the  flow  gradients  are  smail  and  the  fluid  flow 
velocity  is  not  too  supersonic.  The  presence  of  counterpropagating  axial  supersonic  flow  reduces  the 
amount  of  filamentation  growth,  so  the  results  presented  here  may  overestimate  filamentation. 

The  continuity  and  momentum  equations  for  the  quasi-neutral  mo-fluid  plasma  .n  ' h •*  pre«<*ne-' 

of  the  laser  ponderomotive  force  are  linearized,  and  combined  to  give  a  Iriven  ion-  ac  >ustic  wave 

JO 

equation  for  the  electron  density”  : 
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where  a=Ze  end  C^  =  (ZTe+T.)/m.,  and  a  phenomenological  damping  term  V  has  been 

included.  This  equation  describes  the  plasma  responding  as  an  ion-acoustic  wave  driven  by 
temperature  (pressure)  gradients  (first  term  on  the  right)  and  the  ponderomotive  force  (second  term). 
The  variable  ln(n  )  preserves  the  correct  nonlinear  isothermal  steady  state  behavior,  and  ensures 
positivity  of  the  density  in  the  transient  regimes.  The  term  ~V ^ln(ne)V  i  is  ignored,  as  it  is 
second  order  in  the  perturbation.  For  all  of  the  cases  considered  in  this  paper,  the  condition 
Jn^/n  <<1  is  valid;  typical  values  of  5n J n  in  the  time-dependent  calculations  are  less  than  ~5%. 
(For  f§\  laser  light,  hot-spots  shift  randomly  on  the  order  of  the  coherence  time  (~psec),  which  is 
faster  than  than  the  density  can  respond.  Intense  LSI  filaments  are  not  precluded,  though,  since  even 
shallow  density  channels  can  produce  substantial  refraction  over  long  propagation  distances.  For  long 
pulse  non-ISI  lasers,  quasi-steady  state  density  equilibration  can  occur,  and  a  different,  nonlinear 
steady-state  formulation  for  the  electron  density  is  used.  This  algorithm  will  be  described  later.) 

The  ion  temperature  is  assumed  to  be  constant  in  time,  since  the  ion-electron  energy 
equilibration  time  is  typically  on  the  order  of  nanoseconds  in  these  plasmas,  and  this  is  much  longer 
than  time-scales  of  interest.  The  background  ion  and  electron  temperatures  are  assumed  to  be  equal. 
(The  model  and  results  are  insensitive  to  the  ratio  of  background  electron-to-ion  temperatures,  and  in 
any  case  the  results  can  be  easily  renormalised  to  account  for  different  ratios.)  The  perturbed 
electron  temperature  is  then  found  by  solving  the  relevant  energy  balance  equation: 
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Q  is  the  electron  thermal  heat  flux  -C^V S  is  the  Joule  heating  source  given  by  li-^l  /C^TJ 
is  the  electron  thermal  conductivity  (including  any  flux-limiting  effects),  and  ^(n^TJ  tlle  i“verse 

bremsstrahlung  absorption  coefficient.  Compression  effects  and  electron-ion  energy  coupling  are  ignored 
since  they  are  much  smaller  than  the  terms’  included  in  (2.5)  for  the  cases  studied  here. 

The  equations  (2.4)  and  (2.5)  can  be  put  into  a  form  that  depend  on  dimensionless  quantities 
describing  the  magnitudes  of  relevant  physical  phenomena.  We  define: 

ponderomotive  pressure 


plasma  thermal  pressure 
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thermal  conduction  transit  time  across  laser  wavelength  dimensions 
ion-acoustic  transit  time  across  laser  wavelength  dimensions 
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inverse  bremmstrahlung  heating  rate 


Tt2  = 


thermal  conduction  cooling  rate  across  laser  wavelength  dimensions 
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where  A  and  2  are  the  plasma  ion  mass  and  charge,  InA  is  the  Coloumb  logarithm,  and 
^(Z)=(Z-t-.24)/(l+.24Z).^  The  subscript  ’o’  again  refers  to  evaluation  of  the  variable  at  an 
unperturbed  or  initial  value. 

is  the  scaling  constant  of  the  steady-state  temperature  perturbation  in  the  plasma,  which 
is  the  active  force  in  the  thermal  filamentation  mechanism.  The  sensitive  dependence  upon  plasma 
temperature  is  due  to  the  temperature  dependencies  of  the  electron-ion  collision  frequency  (~T^  ^ “) 
and  electron  thermal  conductivity  (~T^'  ).  As  the  temperature  rises,  both  the  decreasing  collision 
frequency  and  increasing  thermal  conduction  act  together  to  smooth  temperature  gradients  and  quench 
the  mechanism.  For  this  reason  thermal  filamentation  is  more  important  in  cooler  dense  plasmas, 
such  as  those  created  by  short  wavelength  (XQ<0.5/lm)  irradiation.  The  ponderomotive^force,  on  the 
other  hand,  is  more  important  in  long  wavelength  laser  irradiation  of  plasmas  (7  ~IoXo/T  )  because 
of  its  explicit  dependency  on  laser  wavelength.  r 

Another  important  difference  between  these  two  types  of  filamentation  is  the  mechanism  by 
which  they  couple  the  laser  light  to  the  plasma.  The  ponderomotive  force  is  almost  instantaneously 
felt  by  the  plasma  (on  times  of  order  w’  ),  and  is  stronger  for  hot  spots  with  shorter  scalelengths. 
In  contrast,  the  thermal  forces  in  the  plasma  require  the  establishment  of  temperature  gradients.  The 
temperature  distribution  is  created  on  a  characteristic  time  scale  that  is  longer  for  larger  scalelengths 
(T~neL  /*e),  so  the  force  is  transmitted  more  slowly  for  large  scale  filaments.  At  shorter 
scalelengths,  the  force  is  suppressed  by  the  diffusive  thermal  conduction  smoothing.  Thus,  thermal 
forces  are  greatest  at  large  scalelengths  and  long  times,  as  opposed  to  ponderomotive  forces  which  are 
greatest  at  short  scalelengths  at  all  times.  These  characteristics  are  the  basis  of  the  qualitative 
differences  in  the  two  filamentation  mechanisms. 

Using  the  quantities  (2.0),  the  plasma  response  equations  (2.4)  and  (2.5)  can  be  re-written  as: 
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Again,  the  spatial  coordinate  x  is  normalized  to  the  laser  wavelength  X  and  the  time  is  normalized 
to  X  /C  ,  the  ion-acoustic  transit  time  across  a  laser  wavelength.  Tlie  overbar  indicates  that  the 
varia&le  is  normalized  with  respect  to  its  initial  or  unperturbed  value,  e.g.,  n  =  n  (x,z,t)/n  (z)  and 

T  =  i(^*,t)/ro(*=o). 

In  the  computer  code,  the  plasma  variables  n^x.s.t),  3ne(x,z,t)3t,  and  Te(x,z,t)  are  defined  on 
the  same  {x,z}  mesh  as  the  laser  electric  field.  The  numerical  algorithm  for  the  hydrodynamics  of 
eqn.  (2.7)  uses  a  combination  of  FFTs  and  an  analytic  solution.  Assuming  relatively  small  variations 
in  the  sound  speed  in  x,  the  driven  ion-acoustic  wave  equation  can  be  Fourier  transformed  to  yield: 
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where  q  is  the  ratio  of  the  imaginary  to  real  frequency  of  the  ion  acoustic  wave,  j  ('  !:i  :h-* 

analysis  presented  here,  q  is  generally  taken  to  be  1/2,  in  accordance  with  the  iigtuiicoiu  m  ivi 


damping  for  Te~T^).  For  the  purposes  of  the  calculation,  the  driving  terms  are  assumed  to  be 
approximately  constant  over  the  time  step  used  in  the  code  (on  the  order  of  a  picosecond).  The 
solution  to  (2.9)  for  a  source  that  is  constant  from  t  to  t  +7"  is  then  found  analytically: 
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where  (1-q  /4)  (  and  G=  {7pi*  ^(k,1t,+T‘/2)-rT(kx,t+7‘/2)/(l-l-l/Z)}.  This  result  is  inverse 
Fourier  transformed  numerically  to  yield  the  plasma  density  at  the  advanced  time  t  -f -T. 

The  plasma  energy  balance  eqn.  (2.8)  is  solved  with  a  three-point  predictor-corrector  method 
and  is  subject  to  the  same  periodic  boundary  conditions  imposed  by  the  Fourier  treatment  of  the 
laser  light  and  plasma  density  calculations.  This  can  cause  problems,  since  the  net  energy  deposited 
by  inverse  bremmstrahlung  at  a  given  axial  position  cannot  be  lost  either  by  heat  flow  out  of  the 
system  or  conversion  to  axial  plasma  kinetic  energy.  Under  these  circumstances,  the  plasma  at  any 
given  axial  position  will  continually  gain  energy  (temperature)  as  the  interaction  progresses,  in 
contrast  to  the  real  system  which  will  reach  a  qnasi-steady  equilibrium  with  laser  energy  converted  to 
plasma  blowoff  and  target  acceleration.  To  avoid  this  unphysical  behavior,  the  spatially  averaged 
energy  gain  (4T2/dx>7T2(x-,t*)f  (x^t’j/X^^)  is  subtracted  from  the  source  term  of  eqn.  (2.3)  to  keep 
the  mean  temperature  constant.  This  is  equivalent  to  assuming  that  the  energy  losses  (to  axiai 
conduction  or  flow)  are  independent  of  the  coordinate  x. 

A  harmonic-mean  flux-limited  formalism  is  used  for  the  heat  flux: 
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is  the  Spitzer  electron  thermal  conduction  coefficient  of  the  plasma  and  f  is  a  phenomenological 
flux-limiter.  Anomalous  flux  limiting  in  laser  plasmas  has  been  implied  by  indirect  experimental 
measurements  of  axial  heat  flow  in  the  region  between  the  critical  density  and  the  abiac-.on  surface. 
However,  the  magnitude  of  flux  limiting  in  transverse  heat  transport  in  the  underdense  region  is 
speculative  at  best.  The  semi-classical  value  f=0.1  is  used  here  simply  to  ensure  that  the  heat  flux 
stays  within  physically  allowed  bounds'  .  In  most  of  the  calculations  presented  here  (and  in  ail  of 
the  ISI  calculations)  the  heat  flux  is  much  smaller  than  this  limit,  so  the  results  are  insensitive  to 
the  flux-limiter. 

For  time  independent  problems  (i.e.,  loser-plasma  interactions  with  non-ISI  lasers  and  puise 
lengths  long  compared  to  hydro  times),  the  calculation  of  the  plasma  density  is  simplified.  The 
solution  of  the  steady-state  momentum  balance  equation  for  time-independent  sources  f  ^(x,z)  and 
Te(x,t)  is 


ne(X’2)  = 


Cg(X,Z) 


•«s 


Cs  (X*,t)  Vx„*  f(x’,z)) 


(2.12) 


where  C’  is  a^constant  of  integration  given  by  the  definition  of  average  density:  C'  = 

n  X _ ’  rdx’C  "(x*)»xp(-7  fdx’C  ”7  ,9  9).  The  temperature  T  (x.z)  :s  found  by  iterative  solution 
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ot  eqn.  (2.3)  with  arbitrary  time  steps,  performed  until  the  solution  converges.  These  separate 

solutions  for  T^fx.i)  and  n^(x,z)  are  then  iterated  alternately  to  obtain  a  convergent  steady-state 
solution  for  both. 


HI.  Analytic  Formulation  of  Filamentation 

Before  proceeding  with,  the  computational  analysis  of  filamentation,  we  present  an  analytic 
treatment  of  filamentation  that  is  based  upon  the  perturbation  solution  of  the  complete  Maxwell’s 
wave  equation.  The  formalism  is  developed  in  general  terms,  then  applied  to  the  ponderomotive  and 
thermal  filamentation  mechanisms  separately.  The  general  formalism  is  also  heuristically  extended  to 
include  spatially  and  temporally  incoherent  ISI  light.  Results  for  both  ponderomotive  and  thermal 
filamentation  of  ISI  light  are  also  presented.  In  section  IV,  we  will  then  compare  results  of  the 
computational  and  analytic  treatments  of  filamentation. 

A.  Basic  theory  of  Filamentation 

The  formalism^  used  here  to  describe  the  scaling  of  thresholds  and  growth  rates  of  filamentation 
is  well  established'”  .  Perturbations  of  plane  waves  in  homogeneous  plasmas  are  analyzed  using  the 
exact  Maxwell’s  wave  equation  for  the  electric  field  and  an  expansion  of  the  dielectric  constant  in 
terms  of  the  perturbed  laser  intensity,  £(<(E+5E)^>)-6  +e’E*5E.  Choosing  E=  E  x*^  \-k  z), 
and  as^suimng  the  perturbed  field  5E  has  the  form  5^Qexpji(Uot-k^z)  +  k  z  +  i(k^x  — k  y)j°  where 
k  <<k  +k“,  one  finds  a  general  dispersion  relation  for  both  ponderomotive  an!  thermal  mechanisms: 
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where  £€  =  9£(}f^Jk^JIo-f<5l)/6[<5l/IQ];  5£  =  (n/nc)7  for  the  ponderomotive  mechanisni  and 
<j£=(n/n JT^/I^x^kyl  the  thermal  mechanism.  fFor  notations!  convenience,  we  define  k  as  a 
wavenumber  normalized  by  the  laser  vacuum  wavenumber;  i.e.,  k|=k^c/Wo).  We  have  assumed  5E 
parallel  to  E,  as  this  is  the  fastest  growing  configuration.  For  the  one  dimensional  case  we  let 
ky=0,  kx=k|,  and  assume  that  the  two  dimensional  case  is  confined  to  k^  =  k  =k|/V-;  then  the 
general  expression  for  the  instability  threshold  in  N  transverse  dimensions  (where  1^=1  or  2)  is: 
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For  any  interaction  strength  (5e ) ,  there  will  be  some  range  of  unstable  perturbation  wavelengths  of 
the  incident  light  intensity.  There  remains  the  questions  as  to  whether  thi3  range  is  relevant 
(contained  in  the  interaction  region),  and  if  so,  whether  the  perturbations  havp  room  to  grow  within 
the  propagation  region.  The  spatial  growth  rate  has  a  maximum  at  a  value  k™a"X  determined  by  the 
a  root  of  the  equation: 
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where  ff£’  is  differentiation  of  <5f  with  respect  to  (k|)  For  the  ponderomotive  mechanism  5e'-0, 
and  for  the  thermal  mechanism  <?£’=  (n/n c)7'r*>-  Using  these  relations  in  (3.3)  yields  the  fastest 
growing  modes  and  their  associated  growth  rates  for  the  thermal  and  ponderomotive  cases.  The 
results  are  shown  in  Table  I;  they  will  be  compared  to  numerical  solutions  later. 

Although  the  fastest  growing  thermal  mode  has  an  infinite  wavelength,  the  growth  rate  is 
approximately  constant  for  k^<<l,  which  is  generally  the  region  of  interest.  Thus,  for  thcrrqai 
filamentatipn,  most  modes  grow  at  the  same  rate.  In  contrast,  the  ponderomotive  grqwch  rate  is  ~k  i 
for  small  k^.  If  the  ponderomotive  force  is  strong  enough,  there  is  some  value  of  k|  at  which  the 
ponderomotive  and  thermal  growth  rates  are  equal.  This  point  occurs  for  filaments  of  size 
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A  i /\„=(7p/7'po)  *  Ponderomotive  filamentation  is  dominant  for  filaments  smaller  than  this  si: 

wnile  thermal  liiamentation  is  more  important  for  the  larger  filaments.  This  suggests  a  '■•"•ur 


wnue  tnermai  .uamentation  is  more  important  ;or  the  larger  ::iaments.  inis  suggests  a  '.-"nir: 
whereby  the  thermal  and  ponderomotive  mechanisms  can  coupic:  large  wavelength  modes  can  :  egm 
focus  to  due  thermal  liiamentation.  then  become  iominated  by  pomierotnet :ve  filamentation  vh-n  ‘h-’ 
filament  becomes  small  -notivh. 


B.  ISI  and  filamentation. 

The  ISI  optical  smoothing  technique  leads  to  rapid,  random  fluctuations  in  the  laser  intensity. 
If  we  average  these  fluctuations  over  a  time  t  ,  then  the  fluctuations  compared  to  the  mean  are 
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proportional  to  (t  /t  )  '  ,  where  is  the  coherence  time  of  the  laser.  Since  the  characteristic 
plasma  hydrodynamic  response  time  is  much  longer  than  the  typical  laser  coherence  time,  the  plasma 
dielectric  constant  will  respond  to  the  time-averaged  intensity.  As  a  simple  way  to  account  for  this, 
we  will  substitute  this  time-averaged  intensity  perturbation  (IQ[tc/t  ]  )  f°r  *ke  background 

intensity  (Iq)  in  the  formulas  for  the  dielectric  response,  using  a  suitable  value  for  t&  .  After  this 
substitution,  the  filamentation  formulas  given  in  the  last  section  will  be  re-derived. 

This  analysis  ignores  some  effects  which  may  contribute  to  the  suppression  or  enhancement  of 
filamentation.  For  instance,  the  stochastic-like  fluctuations  in  the  density  will  increase  the  light 
scattering  and  counteract  filamentation.  On  the  other  hand,  the  time  average  of  the  intensity 
fluctuation  is  treated  linearly,  althongh  the  interaction  is  itself  nonlinear;  this  underestimates  the 
filamentation.  It  is  implicitly  assumed  that  corrections  due  to  these  effects  are  small. 

The  averaging  time  t  is  taken  to  be  the  characteristic  time  for  the  filamentation  mechanism 
to  change  the  dielectric  constant  appreciably  over  the  transverse  dimensions  of  the  filament,  Xi .  For 
the  ponderomotive  mechanism  this  averaging  time  is  the  ion-acoustic  transit  time  of  the  filament. 
X|/C  .  For  the  thermal  mechanism,  the  averaging  time  is  the  larger  of  the  ion-acoustic  transit  time 
and  tile  thermal  conduction  transit  time,  ne^j/^e-  The  ion-acoustic  transit  time  is  larger  when  the 
ratio  T-n^lAo  (3ee  e<ln  (2-6b|)  **  l*33  than  one.  This  ratio  is  largest  for  cooler,  short  wavelength 
plasmas  17,^,  ~l/(T^Xo)}  and  large  filaments.  In  a  worst  case  (Xo=0.25/im  in  a  CH  plasma  with 
T  =  lkeV,  n/n  =0.25)  the  ion-acoustic  time  is  the  dominant  averaging  time  for  X^£350/Im.  Since 
this  worst-case  value  of  350/Jm  is  larger  than  almost  all  filaments  studied  here,  we  will  use  the 
acoustic  transit  time  as  the  averaging  time  for  both  thermal  and  ponderomotive  filamentation. 

The  filamentation  analysis  presented  in  section  III. A  is  now  repeated,  except  that  the 
background  dielectric  constant  depends  upon  the  time  averaged  intensity  perturbation  level,  instead  of 
the  plane  wave  intensity.  Thus,  we  substitute  I0(tc/(;av -)1'  ”  (with  tav^  =  X ^/CJ,  for  IQ  in  the 
expressions  for  the  dielectric  change,  ffc,  given  in  section  III. A.  For  the  ponderomotive  mechanism 
this  substitution  gives: 


5e  =  a  7  r1/2[k,]1/2 

p  n  'p  c  L 


fie*  =  -  7  -  7  r1/2[k.]5/2, 

p  4  n  'p  c  L  1J 


where  we  have  defined  a  normalized  coherence  time  as: 


r  .  !c^s  \illim  T  (keV)l 1/2  tc~— - 

Tc  "  X  “  -31  l  A  VkeV'J  XJ/;m)  • 


For  the  thermal  mechanism, 


x  «  n  1/2 r7  ,-3/2 
55 1  = ;  ’t2tc 


<y€'  _  2  n  r1/2f]c  l1/2 
0€t  4  a  7T2c  L*1J 


Using  these  equations  in  eqns  (3.1)-(3.3)  gives  the  relations  for  ISI  filamentation  found  in  Table  II. 
These  results  will  be  compared  to  results  of  the  numerical  simulations  in  the  next  section. 
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IV.  Analysis  and  Discussion 

In  this  section,  we  use  the  code  described  in  section  II  to  simulate  filamentation  behavior  of 
different  laser  beam  profiles.  We  proceed  in  a  carefully  structured  manner,  beginning  with  simple, 
well-characterised  problems  having  known  solutions.  Complexity  is  added  to  the  interaction  a  step  at 
a  time,  gradually  building  up  to  simulations  of  realistic  laser  beams  propagating  in  laboratory-type 
laser-plasmas.  As  we  proceed,  unique  characteristics  of  the  newly  added  features  are  identified  and 
discussed.  The  analytic  results  of  section  III  are  also  compared  to  these  simulations,  and  the 
limitations  of  the  analysis  are  addressed. 

We  begin  with  the  simplest  case  of  filamentation:  the  propagation  of  a  Gaussian  laser  beam  in 
a  homogeneous,  nonabsorbing,  medium.  Analytic  solutions  of  this  model  for  ponderomotive  and 
thermal  filamentation  are  developed  in  the  appendix,  and  are  compared  to  the  numerical  solutions. 
This  comparison  exposes  the  limitations  of  analytic  methods  for  even  the  simplest  cases  of 
filamentation,  and  illustrates  the  physics  that  distinguishes  the  two  types  of  filamentation.  Next,  the 
beam  profile  made  more  complex  while  the  plasma  remains  simple:  we  consider  generic  and  optically 
smoothed  RPS  and  ISI  laser  beams.  These  more  realistic  beams  are  composed  of  randomly  phased 
perturbations  of  different  wavelengths,  and  introduce  the  possibility  of  nonlinear  mode  coupling. 
Finally,  we  add  inhomogeneity  and  absorption  to  the  plasma  model,  and  consider  laboratory  plasmas 
relevant  to  ICF  experiments. 

Throughout  this  analysis,  we  will  use  two  basic  parameters  to  measure  filamentation  effects:  the 
filament  focal  length  and  the  focal  intensity  maximum.  These  quantities  are  useful  filamentation 
measures  for  two  reasons:  first,  these  parameters  can  be  directly  compared  to  the  analytic  theory. 
Secondly,  the  focal  length  gives  the  minimum  size  of  a  plasma  in  which  filamentation  effects  may  be 
observed,  while  the  maximum  intensity  quantifies  the  impact  of  filamentation  on  other  nonlinear 


processes. 

The  filament  focal  length  and  intensity  maximum  are  not  easily  defined  except  in  the  simplest 
of  cases.  In  these  simple  cases,  the  propagation  of  a  single-peaked  incident  intensity  profile,  the  first 
intensity  maximum  in  the  propagation  direction  corresponds  to  the  focus  of  the  filament.  In  more 
complex  beam  profiles,  however,  there  may  be  (and  usually  are)  many  intensity  maxima  along  the 
propagation  distance,  of  varying  degrees  of  magnitude.  The  usual  practice  here  is  to  identify  the 
focus  at  the  first  intensity  maximum  encountered  along  the  propagation  direction;  further,  it  is 
supposed  that  this  first  maximum  is  due  to  focusing  of  the  fastest  growing  mode.  In  some  cases, 
however,  Intensity  maximums  further  along  in  the  propagation  may  be  considerably  more  intense  than 
the  first.  In  this  case,  more  than  one  focal  length  or  intensity  maximum  can  be  defined. 

In  the  following  discussion,  we  make  the  following  definitions  for  notational  convenience:  the 
maximum  value  of  a  distribution  I(x,*,t)  over  all  values  of  the  variable  x  is  denoted  os  MAX{I}  , 
while  the  first  maximum  of  the  distribution  in  the  direction  of  the  variable  i  (i.e.,  where  dl/dx=0 
and  d  I/d*  <0  for  the  smallest  value  of  x)  is  denoted  MAX^{I}j.  Also,  the  filamentation  focal 
length  is  sometimes  abbreviated  as  1- 

A.  Filamentation  in  Homogeneous,  Nonabsorbing  Plasmas:  Gaussian  filaments 

The  accuracy  of  theoretical  approximations  and  predictions  is  evaluated  using  the  laser-plasma 
propagation  code  previously  described.  For  the  purposes  of  the  calculation  and  comparison,  we  begin 
by  using  a  nonabsorbing  and  homogeneous  plasmv  Nonabsorbing  means  that  the  laser  energy  is  not 
depleted  as  It  propagates  through  the  plasma;  however,  the  laser  is  allowed  to  heat  the  plasma  to 
produce  the  thermal  and  pressure  gradients  needed  for  thermal  filamentation.  Homogeneity  refers  only 
to  the  background  plasma;  again,  the  laser  is  allowed  to  produce  the  inhomogeneity  needed  for 
filamentation.  Although  unrealistic  in  some  ways,  this  plasma  model  provides  a  good  test-bed  for 
basic  filamentation  phenomena;  it  shows  the  qualitative  filamentation  behavior,  and  allows  us  to 
compare  the  calculations  to  the  non-absorbing,  homogenous  plasma  theory  of  section  III.  The  first 
order  effects  of  absorption  and  inhomogeneity  only  alter  the  quantitative  behavior  of  filamentation. 
Absorption  counteracts  filamentation  to  first  order  by  decreasing  the  growth  rate  by  the  amount  -K. 
inhomogeneity,  on  the  other  hand,  causes  local  variations  in  the  strength  of  the  focusing  forces. 
Inclusion  of  these  factors  is  considered  in  section  F/.B,  which  addresses  laboratory  plasmas. 

We  begin  by  considering  the  focusing  of  a  Gaussian  filament  in  the  steady  state  limit.  The 
numerical  results  can  then  be  compared  directly  to  approximate  solutions  of  the  steady-state  aoniin-ar 
propagation  equation  (-.3).  The  derivation  and  results  of  this  Gaussian  model  for  the  ponderomotive 
and  thermal  mechanisms  are  given  :n  the  appendix  to  this  paper:  the  -caiing  and  magnitude  of  the 


focal  lengths  found  there  are  the  same  or  similar  to  the  results  given  in  Table  I.  The  Gaussian 
model  also  predicts  peak  intensities  reached  during  self-focusing.  We  compare  the  calculations  to  the 
results  of  this  model  for  the  cases  of  ponderomotive  and  thermal  filamentation  separately. 

As  a  first  example,  we  calculate  the  intensity  as  a  function  of  the  distance  of  propagation  into 
the  plasma  for  a  case  where  only  the  ponderomotive  force  is  active  (fig.  1).  The  initial  radius  of  the 
filament  is  aQ=20Xo  and  the  interaction  strength  is  7  =5X10  ®  (this  is  40  times  the  threshold  value). 
The  model  and  calculation  results  are  comparable:  the  focal  length  (350\o)  lies  within  ~10%  of  the 
predicted  result  (315X  ),  implying  that  the  Gaussian  model  is  fairly  accurate  in  determining  focal 
.lengths;  this  agrees  witn  previous  findings^'*^.  However,  the  predicted  focal  intensity  is  much  larger 
than  we  find  with  the  code.  In  part  this  is  due  to  the  sensitivity  of  the  model’s  peak  intensity  to 
the  focal  spot  radius:  small  changes  in  the  radius  give  rise  to  large  changes  in  the  peak  intensity 


when  the  focal  spot  radius  is  small  (MAX{I}x"l/width).  In  reality,  the  filament  is  not  constrained 
to  stay  Gaussian,  and  the  peak  intensity  is  not  so  sensitive  to  the  filament  radius.  (The  inability  of 
the  Gaussian  model  to  account  for  saturation  effects  in  the  dielectric  response  is  not  important  here, 
as  the  maximum  dielectric  change  in  the  simulation  is  <0.5%). 

After  the  first  focus,  the  filament  usually  behaves  quite  differently  from  the  model;  propagation 
can  be  periodic,  although  it  usually  has  a  more  complex  periodicity  than  predicted.  The  behavior 
appears  to  depend  on  the  power  level  of  the  filament.  Looking  at  the  long  distance  behavior  of  the 
simulation  just  described  (fig.  2b,  40X  threshold  power),  we  observe  the  beam  breaking  up  into  two 
off-axis  filaments  which  focus  twice  independently,  and  then  combining  again  to  form  a  single 
filament.  The  behavior  is  repeated  again,  but  each  time  the  central  single  filament  is  degraded  in 
power;  the  periodicity  is  only  approximate  since  the  original  Gaussian  is  not  reproduced.  At  higher 
powers  (fig.  2c,  450X  threshold),  the  beam  breaks  into  many  filaments,  which  spray  outward  from  the 
region  of  the  first  focus;  in  contrast  a  lower  power  example  (fig.  2a,  10X  threshold)  exhibits  the 
simple  oscillation  predicted  by  the  model,  although  the  oscillation  is  slowly  damped  by  light  escaping 
the  filament.  In  general,  the  number  of  filaments  formed  after  the  first  focus  is  an  increasing 
function  of  the  incident  beam  power1®. 

The  scaling  behavior  of  the  spatial  growth  rate  compares  favorably  with  the  model.  Fig.  3a 
plots  the  focal  distance  as  a  function  of  7p  for  ponderomotively  focused  Gaussian  filaments  with 
initial  1/e  radius  ao=40XQ.  (The  focal  distance  is  defined  as  the  distance  from  the  beginning  of  the 
propagation  in  the  plasma  to  the  point  where  the  intensity  reaches  its  first  maximum.^  A  linear 
least-squares  fit  through  the  data  points  yields  the  empirical  dependence  1^/a.^-J  '  ,  which  is 

comparable  to  the  predicted  dependence  lj./a  -7°’5  (from  the  appendix  and  table  I).  Calculations 
performed  for  filaments  with  a  Gaussian  1/e  intensity  radius  aQ=20X  also  exhibit  this  scaling. 

The  peak  focal  intensities  are  very  weakly  dependent  upon  7  a^/X“,  the  theoretical  controlling 
parjmeter,  an^  are  much  smaller  than  predicted  (fig.  3b^.  Again,  the  approximation 
C(E  E)— £  +  £’E  E  in  the  model  is  not  violated  even  for  the  largest  values  of  7  a"  / X  ” .  This 
indicates  that  deviations  from  the  Gaussian  constant-shape  ansatz  are  the  failure  mode  in  the  model. 
The  deviations  occur  because  the  expansion  £(x)— £(x=0)+.5£^x“  used  in  the  Gaussian  model  is 
substantially  violated  in  the  simulation. 

These  calculations  were  repeated  using  only  the  thermal  filamentation  mechanism;  the  results 
are  shown  in  fig.  4.  The  initial  filament  radius  is  again  a  =40X  ,  and  the  interaction  strength 
parameter  varies  from  7t2~  ^XlO  ®  to  5X10  The  behavior  here  is  not  as  simple  as  in  the 
ponderomotive  case.  At  low  powers  the  focal  length  varies  as  ~r  matching  quite  closely  the 
scaling  of  the  predicted  values,  although  the  magnitude  of  the  focal  length  is  consistently  ~2Q%  larger 
than  the  predicted  value.  Also,  the  maximum  intensities  are  50%-85%  lower  than  predicted,  though 
they  show  the  correct  scaling  with  the  filament  intensity. 

Contrary  to  expectation,  as  the  interaction  strength  increases  to  high  values  (7'j,''>iSdX  10  ),  the 

focal  length  increases  and  the  focal  intensity  decreases.  Closer  examination  of  the  simulation  reveals 
the  cause  of  this  cunous  behavior:  at  high  intensities,  the  temperature  profile  becomes  flat-topped  and 
sharp-sided  (fig.  5).  The  electron  thermal  conductivity  is  a  strong  power  of  the  temperature,  so  the 
hot  region  (high  conductivity)  is  smooth,  and  the  cold  region  (low  conductivity!  has  steep 
temperature  gradients.  The  density  profile  is  shaped  like  the  temperature  profile  (<5n<>~-i5T  ),  so  the 
filament  refract  primarily  at  the  edges.  The  light  refracted  at  this  sharp  edge  is  refracted  into  large 
angles,  and  the  resulting  interference  of  this  light  with  the  main  body  of  the  filament  creates  high 
frequency  intensity  structure.  This  structure  increases  the  net  diffractive  force  of  the  fiiament  ana 


counteracts  the  focusing  tendency  of  the  density  channel;  thus  the  focusing  is  reduced.  The  deliberate 
creation  of  high  frequency  structure  on  the  incident  light  has  been  suggested  as  a  filamentation 
suppression  mechanism0  ,  and  is  also  the  basis  of  the  RPS  illumination  technique.  Here,  however, 
the  structure  is  created  by  the  laser-plasma  interaction  itself. 

The  sharp-sided,  flat-topped  temperature  profiles  are  dependent  on  the  boundary  conditions  of 
the  code.  The  largest  wavelength  temperature  perturbation  is  limited  to  the  size  of  the  numerical 
mesh,  which  is  only  a  few  times  larger  than  the  filament.  Normally,  this  is  of  little  consequence:  if 
there  is  moderate  heating,  the  conductivity  is  approximately  constant  and  the  Fourier  spectra  of  the 
temperature  has  the  same  spectral  range  as  the  source.  When  the  heating  is  strong  however, 
significant  energy  is  transferred  into  both  longer  and  shorter  wavelength  modes  due  to  the 
nonlineanties  in  the  electron  thermal  conductivity  and  collision  frequency.-  This  energy  cannot  be  put 
into  modes  larger  than  the  mesh  size,  so  it  is  forced  into  the  numerically  resolvable  spectrum,  in 
modes  with  shorter  wavelengths.  These  amplified  short  wavelength  modes  can  then  be  further 
enhanced  by  the  temperature  steepening  associated  with  the  nonlinear  conductivity.  Although  these 
particular  results  presented  here  are  due  in  some  part  to  the  computational  constraints,  the  profile 
steepening  effect  is  a  well-known  phenomenon  in  nonlinear  heat  transport.  The  defocusing  effect  for 
strong  thermal  filamentation  exists,  but  the  quantitative  threshold  calculated  here  is  dependent  on  the 
actual  boundary  conditions  of  the  system,  and  is  not  universal.  As  the  size  of  the  system  increases 
(relative  to  the  heated  region),  the  defocusing  threshold  will  tend  to  increase. 

The  results  of  Gaussian  beam  filamentation  show  the  limitations  of  the  analysis  for  ideal 
filaments  in  simple  plasmas.  Ponderomotive  filamentation  focal  length  scaling  with  intensity  agrees 
well  with  theory,  but  the  predictions  of  the  peak  focal  intensity  are  very  inaccurate.  In  addition,  the 
propagation  behavior  is  different  than  predicted;  as  the  interaction  gets  stronger,  more  filaments  are 
formed  as  the  beam  breaks  up  after  the  first  focus.  Calculations  and  predictions  of  thermal 
filamentation  agree  fairly  well  in  scaling  and  magnitude,  but  only  at  lower  powers.  At  higher 
powers,  nonlinear  temperature  variations  cause  focusing  effects  to  weaken.  Aware  of  these  constraints 
and  behaviors,  we  are  prepared  to  examine  filamentation  occurring  with  more  complex  laser 
illumination  profiles. 

B.  Filamentation  in  Homogeneous,  Nonabsorbing  Plasmas:  Complex  Laser  Beams. 

In  this  section  we  investigate  the  behavior  of  realistic  illumination  profiles  in  simple  plasmas. 
Three  different  laser  beam  types  are  considered  here:  the  typical  or  generic  laser  beam,  the  RPS  laser 
beam,  and  the  ISI  laser  beam.  In  part  (i),  these  profiles  are  defined  and  their  characteristic  features 
are  discussed.  Then,  in  parts  (ii)  through  (iv),  each  profile  is  considered  in  turn,  with  the  results  of 
the  filamentation  simulations  presented.  The  qualitative  features  of  these  simulations  are  discussed 
and  the  quantitative  results  are  compared  to  theory. 

i.  Definition  and  Construction  of  Complex  Laser  Beams. 

We  will  first  consider  a  "generic”  laser  intensity  profile,  representing  a  typical,  high-power  laser 
beam.  These  profiles  are  determined  by  many  installation-dependent  (and  time-dependent)  parameters. 
Imperfections  of  the  optical  system  design  or  components,  optical  misalignment,  or  temperature 
fluctuations  present  in  the  optical  components  during  a  specific  shot,  can  cause  unique  aberrations  in 
the  output  beam.  In  gain-saturated  lasers,  the  aberration  structure  is  preferentially  in  the  phase  of 
the  laser  electric  field.  Further  aggravating  the  problem,  the  desired  focal  spot  size  for  large  scale- 
length  laser-plasma  interactions  is  much  larger  than  the  diffraction-limited  spot  size,  and  quasi-near 
field  intensity  distributions  must  be  used.  In  the  quasi-near  field,  even  small  amounts  of  aberration 
show  up  as  significant  structure  in  the  intensity  profile.  Laser  intensity  profiles  are  sometimes 
characterized  by  their  peak-to-valley  intensity  ratio;  !:1  to  10:1  are  not  uncommon  values  for  this 
parameter. 

In  this  study,  a  variety  of  generic  laser  profiles  with  different  root-mean-square  standard 

deviations  are  used  (a^nj=/dx(I(x)-I  v  )“/ Jl(x)~dx).  The  peak-to-average  intensity  difference  is 

times  the  value  of  O  -  so  typica  f'v alues  of  a  are  in  the  range  of  0.25  to  1.0;  0  —0.5  is 

rms’  .  rms  /  rms 

often  used  as  the  representative  value.  In  the  code,  these  profiles  are  constructed  by  aiding 


randomly  phased  electric  field  fluctuations  onto  a  DC  (k  =0)  electric  fieid  term.  Electric  field 
fluctuations  of  all  wavelength}  (excepting  the  DC  term)  have  the  same  amplitude,  and  span  the 
wavenumber  spectrum  (from  lk|l=l/X^^  to  1/2F,  in  increments  of  l,Xr  v,  recail  that  ail  ipaliai 
dimensions  are  scaled  to  \^).  Varying  t'he  k  ^0  amplitude  relative  to  t;ie  DC  amaiitn  ie 
different  values  of  <7  This  construction  gives  two  different  characteristic  am piitud -a  f  r  the 


Fourier  intensity  spectrum  (fig.  8a).  If  the  magnitude  of  the  DC  electric  field  term  is  Aq,  and  the 
magnitude  of  the  finite  wavenumber  E  field  variations  is  a^,  then  the  intensity,  profile  has  Fourier 
amplitudes  that  are  of  order  "*2A<jaj+a?  and  a^,  for  the  wavenumber  ranges  I  It  i  I  ={l/XmaJc;l/2F} 
and  lkjJ={l/2F:I/F}  respectively.  (The  notation  {a:b}  denotes  the  range  of  values  from  a  to  b.) 
For  small  values  of  <J  ,  a,<<A  ,  and  the  larger  wavenumber  range  is  of  considerably  smaller 
amplitude  than  the  shorter  wavenumber  range.  (In  fig.  6a,  with  F/20  optics  and  a  plasma  jnesh 
length  400\  wide,  these  wavenumber  ranges  correspond  to  I  kj.  I  —  {0.0025:0.025}  and  I  kj^  I  = 
{0.025:0.05),  respectively.) 

Next  we  consider  the  random-phase-screen  (RPS)  illumination  smoothing  method.  This  is  a 
potentially  attractive  near-term  method  for  beam  control,  since  it  can  be  used  on  existing  laser 
beams.  The  laser  beam  is  passed  through  a  simple  random  phase  mask  before  it  is  focused  down  on 
the  target.  This  phase  mask  consists  of  a  large  number  of  discrete  areas  which  randomly  apply  a 
phase  shift  between  0  and  2J  radians  to  the  section  of  beam  passing  through  it.  (This  is  slightly 
different  than  the  current  experimental  version  of  RPS,  where  the  randomly  applied  phase  shifts  are 
either  0  or  T.^’**)  At  the  lens  focal  plane,  this  produces  a  high  frequency,  spatially  incoherent 
pattern  with  a  smooth  envelope  determined  by  the  diffraction  pattern  of  the  individual  phase-shifting 
areas.  This  is  similar  to  the  ISI  method,  but  without  the  laser  bandwidth:  the  structure  in  the  RPS 
method  is  stationary  in  time.  In  principle,  if  the  spatial  structure  has  a  high  maximum  wavenumber 
l^mox,  filamentatioa  growth  rate  will  be  suppressed  (similar  to  the  situation  observed  in  section 
rV.A  with  thermal  filamentatioq).  This  principle  can  be  quantified  using  the  relations  in  table  I; 
filamentation  is  stabilized  when  k|£(2ne7  /nc)^'*  (ponderomotive)  or  k^>(2ne7'j>2/nc)^^  (thermal). 
When  there  is  appreciable  energy  in  the  high  wavenumber  modes,  they  diffract  substantially  over 
distances  smaller  than  the  growth  lengths  of  the  unstable  modes.  If  this  small-scale  diffraction  causes 
significant  changes  to  the  structure  or  phase  of  the  larger  unstable  modes,  filamentation  may  be 
suppressed.  Also,  the  presence  of  appreciable  energy  in  these  modes  implies  that  there  is  less  power 
in  the  unstable  modes,  which  contributes  to  stabilization. 

The  RPS  technique  generates  an  intensity  profile  that  contains  both  larger  amplitude  and 
higher  wavenumber  components  than  the  generic  profile.  A  typical  profile  (fig.  6)  has  higher  peak 
intensities  and  higher  spatial  frequencies  than  a  corresponding  generic  profile.  The  wavenumber 
spectrum  of  the  RPS  beam  is  controlled  mainly  by  the  F  number,  of  the  lens:  the  highest 
wavenumber  component  of  the  incident  electric  field  has  a  wavenumber  ki~l/2F.  The  DC  (k^=0) 
electric  field  term  is  of  the  same  order  as  the  finite  wavenumber  (k^0)  electric  field  amplitudes,  so 
the  high  and  low  wavennmber  intensity  ranges  ( I  ki  I  ={l/Xmax:l/2F)  and  I  k  i  I  ={l/2F:l/F), 
respectively)  will  have  amplitudes  of  the  same  order.  Since  RPS  requires  fast  optics  (small  F/ji),  the 
largest  wavenumber  in  these  profiles  is  usually  much  larger  than  in  the  corresponding  generic  profile. 
One  of  the  purposes  here  is  to  bracket  the  acceptable  F/#  range  for  the  RPS  method. 

Finally,  we  also  simulate  the  induced-spatial-incoherence  optical  smoothing  method.  ISI  in  its 
simplest  form  is  produced  by  passing  a  broad-band  laser  beam  through  an  echelon,  or  stepped 
transmitting  plate.  Each  echelon  step,  like  the  random-phase-screen,  produces  a  phase  shift  by 
imposing  a  time  delay  on  the  beam  passing  through.  Unlike  the  RPS  technique,  the  time  delay  of 
each  echelon  step  relative  to  any  other  step  is  longer  than  the  laser  coherence  time,  so  that  the 
spatially  incoherent  structure  produced  at  the  focus  completely  changes  on  the  time  scale  of  the  laser 
coherence  time.  Since  the  coherence  time  (**lpsec)  can  be  made  much  shorter  than  gross  plasma 
hydrodynamic  response  times  (~100ps),  the  plasma  should  respond  hydrodynamically  only  to  the  time 
average  of  the  laser  profile,  which  asymptotically  approaches  a  smooth  envelope  function. 

Of  these  three  laser  profile  types,  only  ISI  is  inherently  a  time-dependent  profile.  For  the 
generic  and  RPS  cases,  the  intensity  profile  is  frosen  in,  and  the  plasma  can  reach  a  quasi-steady 
equilibrium  if  the  laser  pulse  is  long  enough  (multi-nanoseconds).  We  have  performed  time-dependent 
simulations  using  these  stationary  laser  profiles,  and  they  show  close  agreement  with  the  steady  state 
simulations  after  times  on  the  order  of  a  few  X^/C^.  (1°  other  specific  cases,  it  has  been  shown 
that  steady  state  laser  profiles  can  result  in  inherently  time-dependent  behavior^*.  In  our  formalism, 
the  neglect  of  axial  density  coupling  in  eqns.  (2.7)  and  (2.9)  has  ruled  out  the  possibility  of  this 
type  of  nonstationary  behavior,  so  the  issue  is  not  addressed  here.  However,  it  deserves  closer 
investigation  in  the  future.) 


ii-  Filamentation  of  Generic  Laser  Beams 

We  now  consider  the  behavior  of  the  generic  laser  beam  in  the  homogeneous,  aonabsorbing 
plasma  where  only  the  ponderomotive  force  is  operative.  Fig.  7  shows  a  typical  result  (7  =1.3X10'“, 
=.5.  and  F/10  optics):  in  general,  the  hot  spots  behave  like  individual  noninteracting  filaments 
similar  to  the  Gaassian  filament  shown  earlier.  Interaction  between  filaments  tends  to  occur  only  if 
two  hot  spots  happen  to  be  initially  close  (distances  on  the  order  of  their  own  diameter),  or  on  an 
initially  intersecting  path. 

The  behavior  of  the  individual  Fourier  modes  in  this  simulation  can  be  compared  directly  to 
perturbation  theory.  One  might  expect  the  comparison  to  be  adequate  up  to  the  point  where  the 
mode  coupling  is  appreciable.  Fig.  8  shows  the  low  order  Fourier  mode  amplitudes  of  the  intensity 
distribution,  and  compares  them  to  the  results  found  by  applying  the  Factor  exp(k  (k  j ) L )  (f^om  table 
I)  to  the  incident  intensity  distribution.  In  theory,  the  most  unstable  mo<?e  is  at  k  ^  =0.025 
(X.=40\o)  ;  in  the  simulation  this  mode  closely  follows  the  predicted  growth  for  z£400X  ,  after  which 
it  levels  off  and  then  decreases.  Qther  theoretically  unstable  modes  grow  at  rates  different  (and 
generally  less)  than  predicted.  For  kj^0.04,  the  modes  are  supposedly  stable,  yet  significant  energy 
appears  to  be  going  into,  these  modes  after  propagation  distances  of  100-200Xo.  This  shift  of  energy 
from  small  to  large  kj^  is  expected  for  filaments  that  focus  and  contract  in  size. 

Although  the  individual  modes  do  not  behave  as  predicted  due  to  mode  coupling  effects,  the 
focal  length  predictions  are  fairly  accurate.  Fig.  9  shows  the  variation  of  focal  length  and  focal 
intensity  with  7  f°r  simulations  with  F/20  and  F/10  optics  and  <7  =0.5  on  a  400X  wide  mesh; 

the  focus  is  defined  at  the  point  MAX^(MAX{I}^}  .  There  is  close  agreement  between  the  focal 
length  prediction  and  simulation:  the  simulation  results  lie  very  near  to  the  predicted  values  (from 
table  I)  shown  by  the  dotted  line.  (This  is-  partially  fortuitous,  since  the  formulas  in  table  I  are 
independent  of  <7  .  Although  different  values  of  <7  result  in  different  focal  lengths  (e.g.,  smaller 

fljms  gives  larger  ij.],  the  scaling  of  Ij.  with  7  “  similar.)  For  smaller  values  of  7  ,  lr~7"  (table  I) 
since  the  fastest  growing  mode  is  represented  in  the  incident  intensity  profile  (^jnAxl^/F).  For 
7p20-005,  however,  the  fastest  growing  mode  is  limited  by  the  smallest  perturbation  wavelength  in 
the  incident  spectrum:  X|=20X  .  In  this  range  of  7  >  tbe  fastest  growing  mode  is  constant  and  the 
focal  length  scales  as  lf~7p  (table  I).  .In  contrsSt  to  the  focal  length  predictions,  the  intensity 
maximum  predictions  are  very  inaccurate.  Two  regimes  appear  to  have  been  reversed  between  theory 
and  simulation:  at  small  7  ,  we  expect  MAX{I}^z~constant  (since  I  x~7  Xj^  [see  appendixi,  and 
kJl*Jt~7  /  [table  I|),  ana  instead  we  observe  a  rapidly  increasing  For  7^.5*1.005  we 

expect  MLAX{I}_~7p  (table  I,  using  kI[rlB'x~con3tant},  and  we  observe  a  constant  or  slightly  declining 
value  of  MAX{T}x  ^  with  increasing  7  . 

In  summary,  the  focal  length  predictions  seem  to  apply  better  than  expected  for  ponderomotive 
filamentation  of  generic  laser  profiles  composed  of  many  modes.  In  contrast,  the  predictions  of 
intensity  maxima  are  quite  inaccurate  in  the  same  situation.  The  same  tendencies  were  noted  with 
the  Gaussian  beam  simulations  (see  fig.  3). 

With  the  thermal  filamentation  mechanism  acting  alone  on  the  generic  laser  profiles,  we 
observe  a  qualitatively  different  behavior:  the  filaments  tend  to  attract  one  another  over  distances 
much  greater  than  their  own  dimensions  (fig.  10:  7^,,  =  1.2X10®,  F/20,  and  G ^  =0.5).  The  diffusive 
electron  thermal  conduction  produces  temperature  and  density  structures  mucn  larger  than  the  hot 
spots;  the  result  is  the  development  of  large  wavelength  intensity  modes,  composed  of  many  smaller 
scale  filaments.  This  large  scale  *supermode’  undergoes  large  scale  periodic  focusing  and  defocusing 
when  propagating  over  large  distances.  Since  the  most  unstable  mode  Is  at  k|=0,  and  the  growth 
length  decreases  monotonically  as  kj^Q,  the  dimension  of  this  supermode  is  limited  only  by  the  size 
of  the  physical  system.  The  gross  illumination  symmetry  can  be  affected  by  this  filament  clustering 
behavior. 

The  applicability  of  the  perturbation  analysis  to  these  type  of  profiles  is  assessed  by  examining 
a  Fourier  decomposition  of  the  intensity  profile  (jig.  II).  Simulation  and  theory  agree  that  the 
fastest  growing  mode  is  the  smallest  wavenumber  Ik  1  =0.0025)  of  the  system.  The  periodicity  evident 
in  the  lowest  order  mode  matches  up  with  the  gross  periodicity  of  the  distribution  as  seen  in  fig.  10. 
The  theoretical  period  of  this  mode  is  ~5&0\  ,  about  1,3  of  the  measured  periodicity  of  150U\ 
The  growth  rates  for  other  modes  in  the  simulation  also  tend  to  be  smaller  than  predicted,  and  rrt •  1  ■  t «• 
coupling  effects  are  evident  after  propagation  distances  if  \  few  hundred  \  into  the  plasma. 


The  observed  filamentation  parameters  for  thermal  filamentation  of  generic  beams  are  shown  in 
fig.  12  as  a  function  of  7t'»  f°r  beams  with  O  =0.5  and  F/20  optics.  The  theoretical  focal  length 
and  MAXj{MAX{I}x)  values  are  shown  witli  the  dotted  lines.  The  scalings  observed  in  the 
simulation  match  these  fairly  well,  but  the  magnitude  of  the  values  indicate  much  leas  filamentation 
than  is  predicted:  the  focal  lengths  are  "*3X  larger,  and  the  intensity  maximums  are  ~10QX  larger. 
Although  the  intensity  peaks  are  much  smaller  than  predicted,  there  is  a  definite  sensitivity  to 
interaction  strength,  in  contrast  to  the  ponderomotive  mechanism. 

ilL  Filamentation  of  RPS  Laser  Beams 

Next,  we  will  consider  the  effect  of  the  RPS  optical  smoothing  method  on  filamentation.  The 
profile  of  the  RPS  beam,  as  noted  earlier,  differs  from  a  severely  aberrated  generic  laser  beam  mainly 
in  the 'larger  extent  of  the  RPS  wavenumber  spectrum;  this  extent  is  determined  by  the  F/#  of  the 
focusing  lens.  Here,  we  vary  the  F/#  of  the  RPS  profile  and  measure  the  filamentation  intensity 
increase.  These  are  shown  in  fig.  IS  for  four  cases:  (a)  ponderomotive  mechanism  only:  7  =1.2X  iq‘°i 
7'p2=®>  00  ^hermal  mechanism  only:  7  =0,  7<p2  =  1.2XlO'8;  (c)  both  mechanisms:  7  P=1.3X  10'-, 
7T2=*"^*^  •  (d)  no  mechanisms:  7^ =7'J'2=<^‘  There  is  a  clear  trend  towards  sma  l?er  intensity 

maximums  as  the  F/jjl  decreases.  For  optics  SF/5,  there  is  little  increase  in  the  peak  intensity 
compared  to  the  levels  attained  in  free  propagation  (i.e.,  all  filamentation  turned  off).  The  threshold 
perturbation  wavelengths  for  ponderomotive  and  thermal  filamentation  are  about  30X  ,  and  the 
simulations  were  done  on  a  mesh  with  a  transverse  length  of  400X  .  Thus,  as  the  F/yf  decreases 
more  energy  is  put  into  the  filamentjviion-stabilized  wavelength  region,  X|<30Xq.  The  amount  of 
energy  in  the. unstable  modes  (kiik1  =Threshol4  wavenumber)  is  given  approximately  by  the 
expression  (1-k  [F/#|)  ,  so  different  values  of  k1  change  the  magnitude  of  the  filanjsfjtation 

suppression  effect.  Also,  since  ponderomotive  filamentation  tends  to  have  larger  values  of  k 1  than 
thermal  filamentation,  the  RPS  method  should  be  more  effective  in  suppressing  ponderomotive 
filamentation. 

There  are  possible  concerns  with  the  RPS  smoothing  technique:  first,  the  typical  intensity 
maximum  is  still  of  order  10  times  the  average  intensity  value,  even  for  the  fastest  optics.  Another 
concern  is  that  large  laser-plasma  interaction  chambers  and  laser-fusion  reactor  designs  require  slow 
optics  (^F/20)  in  order  to  reduce  damage  to  optics  and  minimise  the  surface  area  taken  up  by  optics 
in  the  interaction  chamber;  fast  optics  can  not  be  used  in  these  applications.  There  is  a  possibility, 
however,  that  the  many  beams  used  in  symmetrically  illuminated  reactor  designs  may  provide  the 
effect  of  fast  optics:  since  any  area  of  the  pellet  will  be  illuminated  by  a  large  number  of  beams 
incident  at  large  relative  angles,  the  small  F/#  (~F/1)  intensity  distribution  that  is  formed  may 
suppress  filamentation. 

iv.  Filamentation  of  ISI  Laser  Beams 

We  now  consider  LSI  illumination  incident  on  the  homogeneous  nonabsorbing  plasma.  Direct 
comparison  with  the  other  methods  is  more  difficult,  since  the  inherent  time  dependence  of  ISI 
irradiation  complicates  the  measurement  of  filamentation.  The  observables,  instead  of  possessing  a 
single  value,  are  now  represented  by  probability  distribution  functions  of  the  independent  variables  2, 
x,  and  t.  As  a  result,  filamentation  parameters  such  as  intensity  peaks  or  focal  distances  can  be 
defined  in  many  different  ways  which  give  different  values;  we  often  use  more  than  one  definition 
when  evaluating  a  parameter.  (In  the  following,  a  time  average  of  a  distribution  I(x,s,t)  is  denoted 
by  the  brackets  <  >,  defined  as:  <I(x,z)>=  /dtT(x,z,t')//dt’.) 

An  example  of  ISI  with  only  the  ponderomotive  force  acting  is  shown  in  fig.  14.  The 
interaction  parameters  are  7^=0.0051,  r.=3.2,  and  ne/n_ri  =0.5;  the  incident  lens  is  F/20,  and  10 
echelon  steps  are  resolved  on  the  200Xq  wide  calculational  mesh.  Fig.  14(c)  shows  the  intensity 
distribution  <I(x,i)>  averaged  over  84  coherence  times.  Little  intensity  magnification  is  seen:  the 
peak  average  intensity  (MAX{<I(x,z)>}^  )  is  only  ~1.5  times  the  incident  peak  average  intensity 
(MAX{ <I(x,i=0) >}x=  1.281^).  One  of  the  most  noticeable  changes  in  the  intensity  is  the 
development  of  high  wavenumbers  in  the  spatial  structure  as  it  propagates  into  the  plasma.  Plots  of 
the  instantaneous  irradiation  at  42  and  34 1^  exhibit  high  intensity  filaments  (peak  intensities  as  large 
as  12X1^)  in  the  plasma.  These  filamentary  structures  move  about  rapidly  in  the  plasma  ar.  i 
produce  a  much  smoother  time-averaged  distribution.  This  is  favorable  for  hydrodynamic  or  other 
long-time  scale  processes  which  will  respond  only  to  the  smoothed  average  intensity  distribution 
However,  many  laser-plasma  parametric  instabilities  (e.g.,  SRS,  SBS,  2'J^,  etc.'  have  growth  ‘.:mes  in 
the  order  of  the  laser  coherence  tune  and  may  respond  to  the  rapidly  shifting  intensity  spike'  7h- 


magnitude  of  the  differences  in  these  two  distributions  is  seen  by  comparing  the  time-averaged 
intensity  maxima  in  the  plasma,  <MAX{I(z,x,t)}^>,  to  the  maxima  of  the  time-averaged  intensity, 
MAX{<I(x,z)>}x  (fig.  15a).  The  averaged  intensity  maxima  show  the  long  term  existence  of  these 
spikes;  the  average  maximum  intensity  is  "101^,  compared  to  the  average  incident  maximum  of  3.81  . 

Another  measure  of  the  importance  of  the  instantaneous  filaments  is  given  by  the  intensity 
probability  distribution  at  a  given  point  in  the  plasma.  The  ('...Tie  averaged)  incident  intensity 
probability  distribution  is  p»(I/IQ)=  exp(-I/Io);  filamentation  effects  cause  enhancement  of  the  higher 
intensity  portion  of  this  distribution  as  the  laser  propagates  into  the  plasma.  When  this 
enhancement  occurs,  it  is  useful  to  know  how  much  energy  resides  in  filaments;  it  could  then  be 
possible  to  estimate  how  much  energy  is  available  to  drive  other  plasma  instabilities.  The  energy 
redistribution  is  found  by  integrating  the  first  moment  of  the  intensity  probability  distribution 
function  to  determine  how  much  energy  is  at  or  above  a  certain  intensity.  We  call  this  the 
integrated  energy  distribution  function  (IEDF),  and  define  it  here  as: 


IEDF(I/Io)  =\  dX  XpI(X)//  J  dX  XPl(X) 


VK 

The  incident  distribution  function  is:  IEDF(I/I  ,z=0)  =  (l+I/I  )exp(-I/IQ).  The  change  in  this 
distribution  is  shown  as  a  function  of  the  propagation  distance  in  fig.  15b,  where  the  high  intensity 
enhancement  due  to  filament  formation  is  evident.  The  total  amount  of  energy  in  the  high  intensity 
region  is  significant:  ~5%  of  the  energy  appears  at  intensities  greater  than  ~9Iq  in  the  bulk  of  the 
plasma.  In  contrast,  at  t—0 ,  ~0.1%  of  the  energy  is  at  intensities  greater  than  9Iq. 

The  density  variations  responsible  for  the  filaments  observed  in  the  instantaneous  intensity 
distributions  are  relatively  small,  MAX{  1  SnJno !  }  ^5%  (fig.  16).  These  shallow  density  channeb 

produce  filaments  because  the  phase  shifts  are  integrated  over  several  hundred  X  propagation 
dbtances.  This  is  particularly  true  in  these  homogeneous  nonabsorbing  plasmas,  since  Ae  coupling  is 
high  over  the  entire  propagation  region.  In  contrast,  the  high  coupling  regions  in  laboratory  plasmas 
are  much  smaller  in  size,  and  occur  after  most  of  the  propagation  (and  most  of  the  absorption)  has 
taken  place;  we  find  in  the  next  section  that  laboratory  plasmas  generally  do  not  give  rise  to  the 
magnitude  of  filamentation  seen  in  these  simple  plasmas. 

The  results  of  this  particular  simulation  can  be  compared  to  the  theory.  Using  Table  II,  we 
find  that  the  fastest  growing  mode  for  7  T1’  2=9.lX  10’3  and  n/n  =0.5  is  X  i  -30X  .  This  mode  b 
close  to  the  minimum  intensity  wavelength  generated  by  the  F/20  optics  (20A  ).  The  characteristic 
growth  length  (k  M  of  thb  fastest  growing  mode  is  X  7180QX  (Table  II),  and  corresponds  to  a  focal 
length  lp320X  fappendix).  There  are  several  ways  measuring  a  time-averaged  focal  length  in  the 
simulation,  of  which  three  are  used  here:  (l.)  the  position  of  highest  time-averaged  maximum 
intensity  (1^  at  MAX{<MAX(I)  >)z);  (2.)  the  position  of  highest  maximum  time-averaged  intensity 
(L  at  MAX{<I>}^  s);  and  (3.)  the  time  -averaged  position  of  maxim  urn  intensity  (1  j.  at 
<MAX(I}x^>).  (I.j’and  (2.)  can  be  found  directly  from  fig.  15(a),  the  third  is  calculated  during 
the  simulation.  These  values  are:  (1.)  200Xq;  (2.)  162. 5Xq;  and  (3.)  242Xq,  respectively.  All  are 
somewhat  less  than  the  predicted  value;  thb  suggests  that  the  averaging  done  by  the  plasma  is 
nonlinear,  as  it  responds  more  to  the  intensity  peaks  in  the  dbtribution  than  the  time  averaged 
values. 

As  a  more  general  test  of  the  ISI  theory,  we  examine  a  larger  range  of  parameters,  and 

calculate  the  scaling  behavior.  Many  ISI  ponderomotive  filamentation  runs  have  been  made  with 

different  values  of  both  interaction  strength  and  coherence  time."'  7  and  T  were  independently 

varied  in  the  ranges  5X10  "*<7  <2X10  and  .05<T  <3.2.  All  simulations  shown  in  fig.  17  were  done 

with  40  echelon  steps,  F/20  optics,  and  n/n^O.5  on  a  mesh  200Xq  across;  the  results  are  not 

sensitive  to  F/#  (for  F^tO)  or  mesh  size.  The  focal  lengths  (using  the  three  methods  of 

determination  mentioned  previously)  are  shown  sj  a  function  of  the  theoretical  scaling  parameter 

The  dotted  line  is  the  focal  length  range  expected  from  theory,  .as  determined  by  the  fastest 

growing  mode  occurring  in  the  incident  radiation.  At  the  lower  powers,  the  fastest  growing  mode 

varies  as  1 "  (table  II).  At  higher  intensities,  the  fastest  growing  mode  in  table  II  has  a 

wavelength  smaller  than  anv  in  the  incident  spectrum  (X  .  =20X  1,  so  X  should  be  the  fastest 

min  ■>'  min  i  ’i-l  ■' 

growing  mode  The  focal  lengths  u  higher  powers  should  cne.n  vary  as  !r-'  7  T '  “.  The  data  :.j 


closer  tf  1/2  variation  throughout;  a  least  squares  fit  using  all  of  the  points  gives 

i  This  scaling  suggests  that  the  dominant  fastest  growing  wavelength  (Xjnix)  is 
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constant. 


The  averaged  maximum  intensities  are  plotted  as  a  function  of  7,  r1/2  in  fig.  17(b).  There  is 
a  distinct  separation  of  these  averaged  values:  the  time  averaged  maximum  values,  <MAX{I}  >,  are 
always  much  larger  than  the  maximum  time  averaged  values,  MAX{<I>}x-  There  seems  to  loe  little 
consistent  variation  in  the  latter  as  the  interaction  strength  is  increased,  whereas  the  values  of 
<MAX{I}x>  increase  as  "(7  T^2]1/3. 

Thermal  filamentation  calculations  were  also  done  using  ISI  illumination  in  these  simple 
plasmas.  An  example  is  shown  in  detail  in  fig.  18;  this  particular  simulation  is  performed  with  F/20 
optics,  Tc=0.53,  n/nc=0.5,  7t2=2**®  ,  and  ,7<p|=S.8X  10”*.  Even  with  incoherent  illumination,  the 
characteristic  signature  of  tnermal  filamentation  is  apparent  in  the  instantaneous  intensity 
distributions:  the  filaments  attract  one  another,  forming  large  modes  composed  of  high  wavenumber 
structure.  Again,  this  is  associated  with  large-scale  density  and  temperature  fluctuations  (Fig.  19). 
As  in  the  ponderomotive  mechanism,  the  time  averaging  smears  out  most  of  the  structure  seen  in  the 
snapshots;  however,  significant  structure  in  both  high  and  low-wavenumber  modes  can  still  be  seen  in 
the  time  average,  especially  in  the  lowest  order  modes. 

Once  again  the  time-averaged  distribution  is  much  smoother  than  the  instantaneous 
distributions  suggest;  the  averaged  maximum  intensity  <MAX{I}  >  is  significantly  larger  than  the 
maximum  averaged  intensity  MAX{<I>}  (fig.  20a).  The  energy  distribution  looks  very  similar  to 
the  ponderomotively  driven  case  (fig.  20b;  cf.  fig.  15b);  the  integrated  energy  distribution  function 
reveals  that  ~5%  of  the  laser  energy  is  at  intensities  greater  than  101  ,  comparable  to  the 
ponderomotive  example. 

In  fig.  21  the  theoretical  predictions  are  compared  with  a  wide  variety  of  runs  at  different 
intensities  and  coherence  times.  7^2  T  wer®  independently  varied  in  the  ranges 

2.SX10‘7<7t,<SX10'6  and  0.5<T?<5.4.  All  of  tfce  simulations  included  the  theoretical  fastest 
growing  mod^  in  the  incident  intensity  spectrum,  so  the  growth  length  should  scale  as 
I  *  '  The  agreement  between  theory  and  calculation  is  better  than  expected  in  both 
scaling  and  magnitude.  Thwe  is  a  large  scatter  in  the  peak  intensity  values  MAX^{MAX{<I>} 
for  small  values  of  because  the  focusing  is  very  mild;  these  intensities  are  only  a  few 

percent  greater  than  the  incident  peak  (fig.  2lb).  Such  peaks  are  probably  due  to  statistical 
scattering  of  the  light,  not  filamentation.  The  expected  separation  between  the  values  of 
MAX{<I>}  and  <MAX{I}  >  is  observed,  and  a  slight  dependence  on  the  interaction  strength 

(~(7T2Ty2jV5)  ^  noted 

These  simulations  of  ISI  in  simple  plasmas  show  that  the  ISI  filamentation  formulas  in  table  II 
are  only  moderately  accurate.  The  focal  length  predictions  fit  the  observed  thermal  filamentation 
behavior  (fig.  21a)  better  than  the  observed  ponderomotive  filamentation  behavior  (fig.  17a).  This 
can  be  explained  by  noting  that  thermal  filaments  experience  longer  time-averaging  than 
ponderomotive  filaments.  The  ponderomotively  unstable  filaments  are  smaller  in  size  (so  the 
averaging  time  Xi/C^  is  smaller)  than  the  thermally  unstable  modes.  Since  they  are  averaged  less, 
ponderomotively  driven  perturbations  have  a  larger  statistical  deviation.  Some  perturbation  modes 
(that  happen  to  be  larger  than  average)  will  grow  faster,  and  produce  shorter  focal  lengths,  than 
predicted.  This  is  the  behavior  seen  with  the  ponderomotive  simulations  in  fig.  17a.  Thus,  it  is 
reasonable  to  expect  our  ISI  filamentation  theory  to  model  thermal  filamentation  more  accurately  than 
ponderomotive  filamentation. 

C.  Filamentation  of  Complex  Laser  Beams  in  Laboratory  Plasmas 

We  conclude  the  analysis  with  a  treatment  of  filamentation  in  more  realistic  (i.e.,  absorbing 
and  inhomogeneous)  plasmas,  in  particular  those  that  are  produced  in  IGF  research  laboratories  or 
foreseen  for  ICF  applications.  We  cannot  attempt  an  exhaustive  description,  as  there  is  a  huge  range 
of  conditions  encountered  in  these  laboratory  plasmas;  instead,  we  concentrate  on  examples  of  plasmas 
that  are  generated  by  moderately  high  power  laser  light  at  wavelengths  of  1.06dm,  0.53dm,  and 
0.25/im.  These  wavelengths  correspond  to  Nd-glass  lasers,  frequency-doubled  Nd-glass  lasers,  and  K_rF 
or  Nd-glass  frequency-quadrupled  lasers,  respectively.  For  the  1.08  and  0.53/lm  cases,  plasma  or-'ldi^s 
were  generated  by  the  NRL  FAST2D  hydrocode  to  simulate  nesx-term  Hat-target  experiments'”'''1  'J,); 
the  two-dimensional  density  and  temperature  distributions  .are  averaged  ;a  the  transverse  direction  to 
produce  the  one-dimensional  density  and  temperature  profiles  used  in  the  code.  The  0  23dm  laser- 
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plasma  profile  was  generated  by  the  NRL  FASTlD  hydrocode0  in  a  simulation  of  a  few-megajoule 
KrF-driven  reactor-siied  pellet  .  The  first  two  plasmas  have  scaleleagths  on  the  order  of  10° X  , 
while  the  last  is  a  plasma  with  a  scalelength  on  the  order  of  10*X  .  These  density  and  temperature 
profiles,  along  with  the  associated  laser  parameters,  are  shown  in  Sg.  22.  In  all  of  the  calculations 
presented  here,  the  laser  light  is  propagated  in  the  region  from  0.01nc  -t  to  Q.5ncrjt;  propagation  to 
higher  density  is  limited  by  the  assumptions  underlying  the  validity  0/  the  parabolic  wave  equation. 
The  filamentation  coupling  is  high  in  the  region  between  0.01  and  0.5ncr.^,  but  the  propagation 
distance  is  small,  and  the  absorption  is  relatively  high  (especially  for  the  shorter  wavelength  plasmas). 
Thus,  little  additional  filamentation  should  occur  in  this  higher  density  region. 

In  each  plasma  there  is  a  range  of  perturbation  wavelengths  over  which  either  the 
ponderomotive  or  thermal  mechanism  dominates.  Thermal  filamentation  dominates  at  larger  scales, 
ponderomotive  at  the  shorter  scales.  As  noted  in  section  III. A,  these  regions  are  delineated  by  the 
characteristic  wavelength  X^/X  =(7p/7,p2^*  •  As  the  lMer  wavelength  decreases  (with  constant 
intensity),  the  relevant  plasma  density  increases  and  the  plasma  temperature  decreases  slightly  (due  to 
the  higher  plasma  heat  capacity).  Thus,  with  smaller  laser  wavelength,  the_  ponderomotive 
contribution  decreases  (7  “X"),  t^e  thermal  contribution  increases  (7'j>->~T^3) ,  and  thermal 
filame_ntation  dominates  over  a  wider  range  of  wavelengths.  In  addition,  the  absorption  rate 
(Av**tteTv”'  )  of  the  plasma  increases  as  the  laser  wavelength  decreases.  The  higher  absorption  raises 
the  filamentation  threshold  and  lowers  the  growth  rate,  since  the  filament  must  now  grow  faster  than 
it  is  absorbed.  We  find  that  absorption  can  effectively  suppress  filamentation  in  shorter  wavelength 
plasmas. 

The  1.06|im  laser-plasma  absorbs  Little  laser  light  in  the  region  O.Oln^  to  O.on  .  ,  and  the 
ponderomotive  mechanism  is  strong  (7p~0.018,  7'j2**“Xl0"<  at  n/n c  ^=0.25^.  The  ponderomotive 
mechanism  dominates  the  thermal  mechsmism  for  X|^270Xo,  which  incfudes  most  of  the  range  of  the 
simulations.  This  plasma  efficiently  filaments  laser  light,  including  incident  light  that  is  very 
uniform.  Fig.  23a  shows  simulation  results  for  generic  laser  profiles  with  different  incident  O __ 
incident  beams  with  perturbation  levels  0  ^£0.05  begin  to  filament.  Since  most  high-power  laser 

beams  have  (7'rms>>0.05,  this  implies  that  fflamentation  is  a  common  event  at  this  laser  wavelength. 
Using  a  random-phase  screen  with  80  phase  shift  sections  and  F/5  optics  does  not  improve  matters 
significantly:  filaments  with  peak  intensities  £251  are  sti"  observed. 

When  ISI  is  applied  at  1.00/im  with  T  =0.25  (t  =Xpsec,  or  Aii//w~0.0 03),  filamentation  is 

appreciably,  but  not  completely,  suppressed  ^compare  Cfig.  24b  to  fig.  23a).  The  instantaneous 

intensity  distributions  (fig.  24a),  the  time  averaged  intensity  majcimums  (<MAX{I)  >,  fig.  24b),  and 

the  integrated  energy  distribution  (fig.  24c)  all  show  evidence  of  instantaneous  filament  formation. 

The  integrated  energy  distribution  reveals  constantly  increasing  levels  of  energy  at  ail  intensity  levels: 

although  the  beam  begins  with  less  than  5X10  of  its  energy  above  10Iq,  at  n/n  =0.5,  more  than 

2%  of  its  energy  is  above  101^.  The  structure  of  these  instantaneous  filaments  (fig.  24a)  reflects  the 

dominance  of  the  ponderomotive  filamentation  mechanism.  The  time  averaged  intensity  distribution 

(averaged  over  250 T ^  (fig.  24d)  is  much  smoother,  since  the  filaments  move  about  and  do  not 

concentrate  in  a  single  area.  Residual  fluctuations  left  on  the  incident  profile  show  little  growth 

compared  to  generic  beams  (MAX{<I>)  <21  in  fig.  24d;  cf.  fig.  23a  where  MAX{<I>}  "10- 
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Filamentation  is  less  dominant  in  the  plasma  created  with  the  Q.SZflm  wavelength  laser. 
Although  a  significant  fraction  of  the  incident  light  gets  to  0.58^-^,  and  the  plasma  is  60%  larger 
than  the  plasma  at  l.OQfim,  the  smaller  laser  wavelength  reduces  the  ponderomotive  force  effect 
(7p-0-0Q57  at  u/n^t=0.25).  Thermal  filamentation  is  stronger  (7-p.,-3X10  '),  and  should  dominate 
ponderomotive  effects  for  filaments  with  X|^80X^.  With  generic  laser  profiles  (fig.  23b),  we  find  that 
tilamentation  begins  to  occur  when  a  moderate  fluctuation  level  (<7_  ^-0. 3)  is  incident,  this 
fluctuation  level  ia  well  within  the  typical  range.  The  RPS  method  with  F  5  optics  suppresses 
filamentatioi'  somewhat,  reducing  filament  intensities  to  about  12  times  the  average  intensity  'square 
symbols,  fig.  23b). 

When  ISI  is  applied  to  the  0.53am  plasma  with  T  =1  (t  _  =  2.  ipsec.  or  Ay  U- 1 .7*  1 0"' ' .  it 
suppresses  filamentation  about  as  weil  as  the  1.00am  case  {fig.  25;  cf.  fig.  24).  Apnafst ;  v 
smaller  interaction  strength  is  oifset  by  the  .ong^r  interaction  region  as  :ri  e  as  u  r”  d  in  .avr 
wavelengths)  in  this  plasma.  Although  thermal  filamentation  is  stronger  m  :hi3  plasma,  the  str'ic.ir' 
of  the  intensity  distributions  suggests  that  ponderomotive  effects  are  stiil  juite  iur.mant. 


Finally,  we  consider  filamentation  in  the  0.25/lm  laser  wavelength  plasma.  At  ne/ncr;t=0-25- 
the  ponderomotive  mechanism  (7  ~8.8Xlo"*)  dominates  the  thermal  filamentation  mechanism 
(7T2~s-7X  10  )  for  filament  sizes  Xi^lSOX^.  Thermal  filamentation  is  actually  less  important  here 
than  in  the  0.53/im  case  because  2  is  smaller  (Z=l  as  opposed  to  Z=3.5)  and  the  temperature  is 
slightly  greater  (due  to  a  higher  intensity  and  longer  pulse  length).  This  example  is  also  different 
from  the  other  two  cases  because  it  has  a  much  longer  plasma  scalelength,  as  would  be  generated  in 
a  direct-drive  laser-fusion  reactor.  The  large  scalelength  presents  a  worse-case  test  for  filamentation 
at  this  wavelength,  since  it  provides  a  longer  gain  path  for  the  unstable  modes.  Counteracting  this 
effect  is  the  higher  absorption  efficiency  of  short-wavelength  laser  plasmas.  The  amount  of 
filamentation  in  these  plasmas  will  be  determined  by  which  one  of  these  two  opposing  effects  is 
dominant. 

Generic  laser  beams  filament  in  this  0.25/tm  plasma  when  the  incident  perturbation  level  is 

greater  than  <J _  ~0.2  (fig.  23c).  The  MAX{I}  vs.  <T  curve  is  similar  to  the  curves  from  the 

rms  .  .  v  ,x,z  rms  . 

longer  laser  wavelength  interactions,  except  that  the  peak  intensities  are  smaller.  These  smaller 
intensities  are  due  primarily  to  the  higher  absorption  rate.  The  RPS  method  with  F/5  optics  does 
not  appear  to  significantly  affect  the  filamentation  tendency  (see  the  square  symbols  in  fig.  23c);  peak 
filament  intensities  £l2IQ  are  still  observed. 

ISI  is  significantly  more  effective  at  suppressing  filamentation  in  this  0.25/tm  plasma  than  it  is 
at  the  longer  laser  wavelengths.  Applying  ISI  with  Tc=l  (tc=0.9psec,  or  &W/W-9X  10"*)  eliminates 
filamentation  over  both  short  and  long  time  averages  (figs.  26a  and  26d).  The  integrated  energy 
distribution  (fig.  26c)  and  MAX{I}^  curves  (fig.  26b)  show  steadily  decreasing  energy  levels  at  all 

intensities  as  the  beam  propagates  into  the  plasma.  In  addition,  the  ratio  MAX{<!>}  - 

MIN(<I>}x)/<I  (z)>  (fig.  28b)  shows  that  the  nonuniformity  level  is  not  increasing  as  the  laser 

propagates.  The  high  absorptivity  of  the  plasma  appears  to  dominate  any  filamentation 
enhancements  due  to  the  longer  plasma  scalelength. 


V.  Conclusions 

We  have  examined  ponderomotive  and  thermal  filamentation  mechanisms  for  Gauaaian,  IS  I , 
RPS,  and  generic  (typical)  laaer  beams  in  laser-produced  plasmas.  Time-dependent  and  steady-state 
laser-plasma  propagation  codes  have  been  constructed  to  simulate  filamentation  under  these  conditions. 
A  standard  theoretical  formulation  of  filamentation  was  presented  and  extended  to  account  for 
incoherent  light,  such  as  that  found  in  ISI  laser  Keams.  The  predictions  of  this  analysis  were  then 
compared  to  the  results  of  the  laser-plasma  propagation  codes.  First,  a  simple  plasma  (homogeneous 
and  non-absorbing)  was  used  to  study  quantitative  aspects  of  filamentation,  and  to  compare  the 
results  to  the  theory.  Then,  simulations  were  done  with  realistic  laboratory  plasmas  to  determine  the 
importance  of  filamentation  in  more  complex  experimental  environments. 

There  is  a  distinctive  behavior  that  differentiates  the  ponderomotive  and  thermal-conduction 
dominated  filamentation  mechanisms.  In  general,  ponderomotively-driven  filaments  interact  locally 
through  interference  effects  of  the  light  waves;  these  filaments  tend  to  be  independent  from  one 
another.  In  thermally-driven  filamentation,  the  high  plasma  conductivity  creates  long-scale  density 
gradients  that  cause  light  filaments  to  attract  one  another  at  large  distances.  This  attraction 
mechanism  decreases  the  spatial  coherence  of  the  beam,  increases  the  width  of  the  perturbation 
wavenumber  spectrum,  and  can  reduce  or  stabilize  further  large-scale  self  focusing.  At  high  powers, 
the  effect  is  enhanced  by  the  nonlinear  behavior  of  the  temperature  profile. 

Simulations  of  Gaussian  laser  beams  show  the  limitations  of  the  theoretical  analysis.  The 
ponderomotive  focal  length  predictions  agree  quite  well  with  the  theoretical  predictions,  but  the 
behavior  of  the  light  in  and  after  the  first  focus  can  differ  markedly  from  the  predictions.  Gaussian 
beams  undergoing  thermal  filamentation  agree  with  theoretical  predictions  only  at  lower  intensities;  at 
higher  intensities  the  focal  length  increases  rather  than  decreases,  and  the  peak  intensities  decrease 
rather  than  increase.  Both  of  these  effects  are  due  to  the  stabilization  effects  of  the  nonlinear 
temperature  profile.  In  both  the  ponderomotive  and  thermal  filamentation  cases,  peak  intensities 
found  in  the  simulations  fall  far  short  of  their  predicted  values. 

There  are  significant  discrepancies  between  the  perturbation  theory  and  the  simulations  for  the 
generic  and  RPS  smoothed  laser  beams:  the  fastest-growing  mode  is  often  different  than  predicted, 
and  the  growth  rates  for  most  longer  wavelength  unstable  modes  are  lower  than  predicted.  The 
supposedly  stable  higher  wavenumber  modes  grow,  apparently  due  to  nonlinear  mode-mode  coupling. 
In  spite  of  these  discrepancies,  comparisons  of  the  results  to  the  theoretical  focal  length  scaling  laws 
show  rough  agreement.  Comparisons  to  the  peak  intensity  scaling  laws  are  again  poor.  The  RPS 
optical  control  technique  is  able  to  suppress  filamentation  under  some  plasma  conditions,  but  requires 
relatively  fast  Incusing  optics  ($F/'5). 

The  ISI  smoothing  technique  is  first  simulated  in  homogenous,  nonabsorbing  plasmas.  Time 
averages  on  the  order  of  a  hundred  coherence  times  show  relatively  smooth  laser  illumination,  buc 
there  can  be  a  simultaneous  increase  in  the  proportion  of  laser  energy  at  the  higher  intensities.  This 
has  important  implications  for  nonlinear  interactions  that  respond  to  the  light  on  times  of  order  of 
the  laser  coherence  time  (~psec).  Scaling  studies  of  the  ISI  focal  lengths  show  agreement  between 
theory  and  calculations.  In  contrast,  the  maximum  intensity  levels  in  the  plasma  are  found  to  be 
relatively  insensitive  to  the  interaction  strength,  in  disagreement  with  our  theoretical  predictions. 

Simulations  using  near-term  laboratory  plasmas  demonstrate  that  filamentation  tends  to  be 
much  stronger  at  longer  laser  wavelengths  (i.e.,  X  =1.06/im  and  0.53/im).  Ponderomotive 
filamentation  is  dominant  at  1.06/lm,  and  is  exacerbated  by  the  relatively  small  underdense-plasma 
absorption.  Generic  beams  of  high  quality  and  RPS-snoothed  beams  both  filament  rapidly 
(producing  filament  intensities  "'25-30  times  the  average).  ISI  also  shows  some  evidence  of  enhanced 
energy  at  higher  intensity  in  the  1.06/im  laser-plasma.  The  strong  filamentation  tendency  exhibited 
by  the  i  Ofi/Jm  wavelength  interaction  (also  at  0.53Um)  underscores  the  importance  of  filamentation  in 
neax-term  laser-plasma  experiments:  the  laser-plasma  interactions  at  these  longer  laser  wavelengths  arc 
probably  dominated  by  filamentation  effects. 


Ia  reactor-sized  plasmas  it  shorter  laser  wavelengths  (X^=0.25/im),  the  increased  absorption 
reduces  filamentatioa  for  ail  types  of  laser  beams.  When  filamentation  occurs,  maximum  intensities 
are  "10-15  times  the  average  incident  intensity.  RPS  smoothed  beams  do  not  reduce  these  maximum 
intensity  levels.  However,  the  ISI  smoothing  technique  or  very  clean  beam  profiles  can 

eliminate  filamentation  in  short  wavelength  laser  plasmas. 
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Appendix 

The  quasi-optic  equation  (2.2)  can  be  solved  using  the  paraxial  approximation,  which  consists 
of  expanding  the  nonlinear  dielectric  constant  in  the  transverse  variable  and  keeping  only  the  first 
order,  quadratic  terms  .  It  is  a  well  known  that  a  beam  that  initially  has  a  Gaussian  intensity 
profile  and  a  constant  or  quadratically  varying  phase  front  remains  Gaussian  as  it  propagates  through 
a  quadratically  varying  dielectric  constant.  We  thus  consider  a  Gaussian  filament  for  consistency. 
The  electric  field  of  the  filament  is  written: 


JKx.7)  =  ?o00 


U(?) 


N/2 


exp{-x2/2a(7)  2+i  (j>Q  (7)  (??)x2)  >  (A.  1) 


where  a(7)  is  the  l/e  radius  of  the  filament  at  distance  TJ,  aQ  is  the  initial  filament  radius,  N  is  the 
number  of  transverse  dimensions,  and  and  ^  are  real.  As  written,  the  field  conserves  energy  as 
a(?7)  varies.  Inserting  (A.l)  in  equation  (2.3) ,  using  the  paraxial  approximation  6(7, x)-  6{tj ,x=0)~ 
(3[e(7lx=0)i/3[x"])x^  (  assume  that  3f/3x  is  real),  separating  real  and  imaginary  parts,  and  equating 
like  powers  of  x,  we  find  the  following  equation  for  the  filament  radius  a (7): 


4x2^-  = 
dr,2 


1  ,  2  3e 

- 3  +  4T  ~ 2 


(A. 2) 


A.l  Ponderomotive  Filamentation 

For  ponderomotive  filamentation,  3e(7,x)/3x2=3(ne(7,x)/nc)/3x2-  a(rj)'2 (ao/a(7/))^ , 

where  7p  is  evaluated  at  x=0.  This  leads  to  the  equation:  ^ 


4r 


2  d2  a.  (v) 


a(7)3 


4t 


(n/a  )7  a 
v  '  c  'p  o 

,  ,1+N 

a(t;) 


(A. 3) 


The  filamentation  threshold  is  determined  when  da(7)/d7=0  and  d  a.(T] ,z—Q) /dT]“  =0  at  7=0.  The 
threshold  filament  radius  is: 


(A.  4) 


This  radius  is  independent  of  the  dimension  N,  and  agrees  with  table  I  if  we 
perturbation  wavelength  of  the  Gaussian  distribution  to  be  X^p(Gaussian) 
(assuming  2n/n</J  <<  N). 

The  solutions  are  dependent  upon  N,  and  are  straightforward  when  n/n 
and  T,=z/e°'  .  For  N=2,  the  solution  of  (A. 2)  is: 


define  the  effective 

=  TVfa  2  4.4a 

o  o 

is  independent  of  7, 


^  (7) 


2 

a 

o 


(l-4,2(./nc)73^) 

4x2e 


(A.  5) 


This  predicts  that  for  filaments  above  threshold,  the  radius  will  go  to  tero  at: 


V'  "J 


The  collapse  of  the  filament  to  tero  radius  at  t=l-,  occurs  in  this  model  because  the  dielectric 
constant  (previous  to  (A. 3))  is  approximated  as  directly  proportional  to  the  intensity.  The  induced  £ 
causes  refractive  forces  that  are  always  larger  than  diffractive  forces,  and  the  filament  collapses  to 
sero  radius  and  infinite  peak  intensity.  In  a  real  plasma,  the  induced  £  saturates  at  £  =  1,  and 
diffraction  eventually  overtakes  the  filamentation  forces  when  the  radius  is  of  the  order  of  X  . 


Although  the  solution  is  technically  invalid  at  s=lrp>  the  collapse  is  rapid  enough  that  the  solution  is 

_ 1**1..*  11  •  1  .  „  i  ••  .  . 


valid  close  to  1fp.  and  ljp  is  a  good  approximation  to  the  focal  distance  in  the  real  case.  This 
value  is  in  agreement  with  table  I  if  we  use  the  previous  definition  of  the  Gaussian  wavelength 
(Xip=TV2ao),  and  define  the  Gaussian  growth  wavelength  X  p=25Tljp. 

When  N=1  the  solution  a(s)  to  (A.3)  is: 


[ca(z)2+2/Ja(z)-l] 


i  r  .  -i 

Tj2  2  *  Sla 


(C  ♦  P2)l/2 


(A.  7) 


_  2xe 


where:  z  =  y^yo  i  C  =  (l-2/Jao)/aQ;  and  p  =  At  (n/n.)7  aQ. 


is  equivalent  to  the  threshold  condition  given  by  (A. 4).  This  displays  a  focal  length  given 


8r*a^(n/n  )rf  £^2 

o  '  cJ  1  p  o 


1/2 
Ta  £  ' 
o  o 


1  /-M-1  'i  _  - - - _ 

£P  (8T2(n/n  )7  -l)3/2  (S?7J1/2 


(A. 8) 


c"  *p  o 


The  approximation  on  the  RHS  of  (A. 8)  is  valid  for  filaments  far  over  threshold:  8T^(n/n  )"]  a^>>l. 
In  this  limit,  lro(N=l)  =  (T/\/8  )ljp(N=2)~  l.lljp(N  =  2) ,  and  matches  table  I  if  we*  define 
V  =4\/?ljp(N=lJ-5.71jp(N=l).  Like  the  perturbation  analysis  in  section  EH,  this  model  predicts  very 
similar  focal  lengths  for  N=1  and  N=2. 

The  minimum  filament  radius,  attained  at  ljp,  is: 


=  [  8r2  -  7  a2  -  1  1 

L  n  'p  o  J 


Thus,  the  maximum  filament  intensity  (MAX{I}x/Iq  =  exhibits  a  linear  dependence  on  the 


incident  filament  intensity. 


A. 2  Thermal  Filamentation 


For  the  thermal  conduction  dominated  filamentation  mechanism,  ;he  dielectric  constant 
n  A  o  1  v 

expansion  gives  Qc{t],x)/Ox“=  9(n  (?7,x)/n  )/9x  -  -2T*(n/n  )7TO(a  /a (tj))1  .  Substituting  this  in 
(A. 2)  yields:  '  ° 


.2  d  afr?)  1  _  4,  ,  .  N  ,  .  1-N 

4t  o  ~  q  8t  (n/a  a(7) 

Q.T]  a(?7) 


(A.  9) 


The  filamentation  threshold  is  then: 

Th  1  f  n  ]"1/4 
aoT  ^  t  la  7T2J 


(A. 10) 


ft 


This  agrees  with  table  I  if  we  define  the  effective  perturbation  wavelength  as  X  j~,=Tv2 the  same 
relation  as  we  have  previously  found.  The  solution  a(7,7-p0)  for  N=1  can  be  found  by  integrating 
eqn  (A.9)  twice.  First  it  is  helpful  to  define  the  transform  variables  u  =  (8T  (n/n^^^a  )  '3a(z) 
and  y  =  (8*4(n/nj7T2ao)2'  s/2*Veo.  Then  (A.9)  can  be  expressed,  in  a  canonical  form: 


dfulzl  _  __1—  _  , 

,  2  “  f  ,3  1 

dy  u(y) 

Integrating  this  twice,  we  find: 

fU^  _ udu 


(A. 11) 


where: 


4y  =  u  du _ 

7  "«•  [2(uo-u)(«-uJ  (u-u_)] 1/2 


*  (8r“  ;  7t2^1/3 

/*  * 


[U8“o] 


This  has  the  solution: 
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+  [2(u  -u  )]1/2  E  -2 -  ;  -2— 1 
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(A. 12) 


where  F($;p)  and  E(j};p)  are  elliptic  integrals  of  the  first  and  second  kind,  defined  as: 

(t  r  2  2  1 -1//2 

F(^;p)  =  dx  (l-pxJ) (1-x*) 

r  o  o  -jl/2 

E(d;p)  =  dx  (l-px“)/(l-x“) 


In  the  limit  u  >>1  (i.e.,  a  >>1),  eqn  (A. 12)  can  be  simplified;  using  the  limiting  forms  E(^;l)  =  ^, 
F(j>;l)=ln(3ec[ain”1^]+tan[3in  1^])J  and  uA=*(2uo)-*' 2  we  find: 


u  1/2 

*  -  -  r la  sot  (''I1-  ^  [«„-»&)] 1/2  ca. i3) 

o  1  w  J  o  J 


The  focus  occurs  where  u(y)=u_j_  (this  is  the  minimum  radius  achieved  by  the  filament),  at 

v  I1/2  K*rlr  w*  fin/J  fru*)il  Unirf.Yi • 


y^(2uo)  '  .  Transforming  back,  we  find  the  focal  length: 


fT  T  (n/n^)  j 


(A.  14) 


This  is  again  independent  of  the  filament  radius,  and  agrees  with  the  sinusoidal  perturbation  result 


(table  I)  if  we  define  the  Gaussian  growth  wavelength  as  X^,p=V2Tlpj,.  The  maximum  intensity, 
which  occurs  at  the  focus,  is: 


(A. 15) 
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a 


References 

A.J.  Schmitt  and  J.H.  Gardner,  J.  Appl.  Phys.  60,  6  (1986). 

S.  Skupaky  and  K.  Lee,  J.  Appl.  Phys.  54,  S671  (1983); 

A.J.  Schmitt,  AppL  Phys.  Lett.  44,  399  (1984). 

Y. R.  Shen,  Prog.  Quant.  Electr.  4,  1  (1975); 

J. H.  Marburg er,  Prog.  Quant.  Electr.  4,  35  (1975);  and  references  therein. 

P.K.  Kaw,  G.  Schmidt,  and  T.  Wilcox,  Phys.  Fluids  143,  1522  (1973). 

C.  Max,  Phys.  Fluids  19,  74  (1976). 

F.W.  Perkins  and  E.J.  Valeo,  Phys.  Rev.  Lett.  32,  1234  (1974); 

V.K.  Tripathi  and  L.A.  Pitale,  J.  Appl.  Phys.  48,  3288  (1977). 

E.  Kane  and  H.  Hora,  Austr.  Journ.  Phys.  34,  385  (1981). 

H.A.  Baldis  and  P.B.  Corkum,  Phys.  Rev.  Lett.  45,  1260  (1980); 

M.J.  Herbst,  J.A.  Stamper,  R.R.  Whitlock,  R.H.  Lehmberg,  and  B.H.  Ripin.Phys.  Rev.  Lett. 
46,  328  (1981); 

Z.  Lin,  Opt.  Conrnmn.  42,  351  (1982); 

K. A.  Tanaka,  B.  Boswell,  R.S.  Craxton,  L.M.  Goldman,  F.  Guglielmi,  W.  Seka,  R.W.  Short, 
and  J.M.  Soures,  Phys.  FL  28,  2910  (1985). 

R.H.  Lehmberg  and  S.P.  Obenschain,  Opt.  Comm.  46,  27  (1983); 

R.  Lehmberg,  A.  Schmitt,  and  S.  Bodner,  Journ.  Appl.  Phys,  in  press. 

Y.  Kato  and  K.  Mima,  Appl.  Phys.  B29,  186  (1982). 

Y.  Kato,  K.  Mima,  N.  Miyanaga,  S.  Arinaga,  Y.  Kitagawa,  M.  Nakatsuka,  and  C.  Yamanaka, 
Phys.  Rev.  Lett.  53,  1057  (1984). 

S. P.  Obenschain,  J.  Grun,  M.  Herbst,  K.  Kearney,  C.  Manka,  E.  McLean,  A.  Mo3tovych,  J. 
Stamper,  R.  Whitlock,  S.  Bodner,  J.  Gardner,  and  R.  Lehmberg,  Phys.  Rev.  Lett.  56,  2807 
(1986). 

Y.  Kitagawa,  T.  Yabe,  Y.  Kato,  and  C.  Yamanaka,  Phys.  Rev.  Lett.  56,  2804  (1986). 

R.H.  Lehmberg  and  J.  Goldhar,  Fusion  Tech.  U,  532  (1987). 

R. S.  Craxton  and  R.L.  McCrory,  Journ.  Appl.  Phys.  56,  108  (1984); 

K.  Estabrook,  W.L.  Kruer,  and  D.S.  Bailey,  Phys.  Fluids  28,  19  (1985). 

S.  Sartang,  R.G.  Evans,  and  W.T.  Toner,  J.  Phys.  D:  Appl.  Phys.  143,  955  (1983). 

J.A.  Stamper,  R.H.  Lehmberg,  A.  Schmitt,  M.J.  Herbst,  F.C.  Young,  J.H.  Gardner,  and  S.P. 
Obenschain,  Phys.  Fluids  28,  2563  (1985). 

A  3D  version  of  the  computer  code  presented  here  has  been  used  to  simulate  ISI  filamentation 
driven  by  the  ponderomotive  force  only.  Preliminary  runs  were  performed  with  1.06/Jm 
wavelength  ISI  light  in  the  plasma  shown  in  fig.  22a.  The  results  indicated  that  although 
some  filaments  grow  to  much  higher  intensities,  they  do  not  contain  much  energy:  the  intensity 
probability  distributions  are  similar  to  the  2D  results  shown  in  section  IV  (fig.  24). 

M.S.  Sodha,  A.K.  Ghatak,  and  V.K.  Tripathi,  Progress  in  Optics.  Vol  XIII  (E.  Wolf,  ed.), 
Amsterdam:  North-Holland,  1976,  and  references  therein. 

A.J.  Schmitt  and  R.S.B.  Ong,  J.Appl.Phys.  54,  3003  (1983). 

M. S.  Sodha,  D.P.  Singh,  and  R.P.  Sharma,  J.  Appl.  Phys.  50,  2678  (1979). 

J. A.  Fleck,  J.R.  Morris,  and  M.D.  Feit,  Appl.Phys.  10,  129  (1976). 

K.  Lee,  D.W.  Forslund,  J.M.  Kindel,  and  E.L.  Lindman,  Phys.  Fluids  20,  51  (1977); 

C. E.  Max,  C.F.  McKee,  and  W.C.  Mead,  Phys.  Fluids  23,  1620  (1980); 

W.M.  Mauheimer,  D.G.  Colombant,  and  J.H.  Gardner,  Phys.  Fluids  25,  1644  (1982). 

N. E.  Andreev,  V.P.  Silin,  and  P.V.  Silin,  Sov.  Phys.  JETP  52,  653  (1980); 

R.W.  Short,  R.  Bingham,  and  E.A.  Williams,  Phys.  Fluids  25,  2302  (1982); 

V.V.  Amelin  and  L.A.  Zelekson,  Sov. J. Plasma  Phys.  10,  625  (1984); 

A. A.  Andreev,  I.G.  Lebo,  and  V.B.  Rozanov,  Sov.  J.  Quantum  Electron.  143,  1412  (1986). 

M.J.  Herbst,  J.  Stamper,  R.  Lehmberg,  R.  Whitlock,  F.  Young,  J.  Grun,  and  B.  Ripin,  NRL 
Memorandum  Report  4983  (1932). 

A.R.  Bell,  R.G.  Evans,  and  D.J.  Nicholas,  Phys.  Rev.  Lett.  46,  243  (1931). 

P.K.  Kaw.G.  Schmidt,  and  T.  Wilcox  Phys.  FI.  1J5,  1522  (L973j. 

V.A.  Petrishev,  Radio.  Quant.  Elec.  14,  1062  (1971). 

D. J.  Nicholas  and  3.G.  Sajjadi,  J.  Phys.  D:  Appl.  Phys.  R),  737  (19861. 


Table  I 


Ponderomotive : 


Steady  State  Filamentatior.  Parameters 
Thermal : 


[Th 


A.  k^  ;  Threshold  Perturbation  Wavenumber : 

©}*  ‘  {£-  *  sH'1 

&  ~  7  o 
n  'p 
c  * 


B.  k  ;  Spatial  Growth  Rate: 

s 


k  p  =  —  k^{2-  7 

&  2fe  il  -  p 


‘-2  _1_ 

£1  '  Ne 


OJ 


->} 


1/2 


kgP  "  kl(2e  n  7p  1) 


1/2 


o  c 


C.  k™^;  Fastest  Growing  Perturbation  Wavenumber: 


E8  ‘  7K} 


1/2 


=  -M2* 

r>rr  1  n 


«T  2J7 


7T2 


1- 


Ne 


k  ~  ( -i  a  7  -  k2] 

gT  l2e  n  7T2  *lj 


o  c 


frmax’l^  1 

r  1  1  r1 

lklp]  *2 

l2n/n  7  +  Ne  J 

’  C'D  O 

kmax  =  0 

“IT 


D  k^11*;  Fastest  Spatial  Growth  Rate: 

o  .-1/2 


max 

1  n  , 

gP  “ 

2  n 

c 

nax„ 

1  n 

'gP  " 

2  a 

J7 


’max  M  n  ) 
gT  l2e  n  7T2i 


1/2 


o  c 


i 


1/2 


mm 


imS 


uare  points  are  calculations  with  aQ/X  — 40, 
il  the  open  circles  are  calculations  with 


(b)  focal  intensity,  I(i=lj),  vs.  ?T2 
Gaussian  beams  with  a  =40\  . 


Fig.  6  Typical  laser  intensity  profiles  for 
(a)  the  generic  laser  beam  (with  F  ''20  optics 
and  ffrmj=0.25);  and  lb)  the  RPS  optical 
technique  (with  F/J  optics).  The  top  figures 
show  the  Fourier  spectra  of  the  eiectnc  fields 
(dotted  lines)  and  the  intensity  (solid  lines); 
the  bottom  figures  show  the  corrcsoonding 
intensity  orofiies  Ifx,i=.1). 


ical  curves  (using  table  1  and  the 
lx);  the  square  symbols  are  the  r< 
optica,  and  the  asterisk  syml 
nils  using  F/tO  optics. 


parameters  are  as  in  fig.  10.  The  dotted  line 
shows  the  theoretical  result  from  table  1  and 
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Fig.  13  Variation  of  MAX{I}x/Iq  vs.  lens 
F  number  shows  decreased  filamentation 
effects  for  fast  optics.  Other  parameters  are 
the  same  as  in  figs.  7  and  10. 


Fig,  14  Ponderomotive  filamentation  of  ISI 
light  in  a  homogeneous  plasma:  contour  ind 
isometric  plots  of  the  intensity  distributions. 

(a)  and  (bj:  instantaneous  intensities  at  42T_ 
and  34T,;  (c)  time-averaged  intensity  vs.  x.i 
over  347"’..  Interaction  parameters  are: 

7  =0.0051,  r  =3.2,  n  ,'n  -,  =  0.5,  and  F'20 

'  n  -  ’  i'  cnt  ’ 
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Fig.  15  ISI  ponderomotive  filamentation  in 
a  homogeneous  plasma:  (a)  Normalized 
maximum  intensities  <MAX{I}x>/Iq  (solid 
line)  and  MAX{<I>}x/I  (dotted  line)  vs. 
propagation  distance,  (b^  Contour  plot  of 
the  integrated  energy  distribution  function, 
IEDF(I/Ia),  or  the  fraction  of  laser  energy 
greater  than  intensity  I/I0i  as  &  function  of 
laser  propagation  distance.  Interaction 
parameters  are  the  same  as  in  Fig.  14. 


■r'-t  :ucu: 


era  ore  the  same  as  in  Fig. 


.r y eat  peak  MAX{^1  AX{<I>)x)t;  (iv.)  the 
oiled  line  allows  the  expected  values  using 
le  results  in  Table  II.  All  simulations  are 
itli  n/n  ..—0.5,  K/20  optics,  on  a  20<)Ao 


Fig.  18  ISI  Thermal  filamentation  in  a 

homogenous  plasma:  contour  and  isometric 

plots  of  intensity  distributions,  (a)  and  (b): 

instantaneous  intensities  at  1257"  and  250T  ; 

,  ,  .  ,  .  c  c1 

(C)  time-averaged  intensity  vs.  x,s  over  250T. 

7t^  =  2X10_6,  r.= 0.53,  7t^=J.3X10'4, 
n  .  n  .  =0.5,  and  F/20  optics. 


Temperature 


Fig.  20  ISI  Thermal  Filamentation  in  a 
homogeneous  plasma:  (a)  <MAX{I}x>/Iq 
(solid  line)  and  MAX{<I>}  /I  (dotted  line) 
vs.  propagation  distance,  z/£q.  (b)  Contour 

plot  of  the  integrated  energy  distribution 
function  IEDF(I/Iq)  (the  fraction  of  energy 
greater  than  a  given  intensity)  as  a  function 
of  the  propagation  distance,  z/X^.  From  the 
simulation  shown  in  Fig.  18. 


max 


Fig.  24  ISI  filamentation  at  1.06/im  loser 
wavelength  in  laboratory  plasma  fig.  2a. 

(a)  Inatantaneoui  intensity  distribution  at 
250T  ;  (b)  time  averaged  intensity  maximums 
<MaSc{I}x>  (solid  line)  and  MAX{<I>)x 
(dotted  line);  (c)  contour  plot  of  ths 
integrated  energy  distribution  function 


(fraction  of  energy  greater  than  a  given 
intensity  I/IQ) ,  IEDF(I/I  jn/nJ;  and  (d)  time 
averaged  distribution  <l(x,s)>. 


Fig.  26  ISI  filamentation  at  0.25iim  laser 
wavelength  in  laboratory  plasma  Tig.  22c. 

(a)  Instantaneous  intensity  distribution  at 
250T  ;  (b)  time  averaged  intensity  maximuma 
<MAX{I)  >  (solid  !  ine)  and  MAX{<I>}x 
(dotted  line),  and  relative  peax-to-vailey 
uniformity  (MAX{<I>)- 
MIN{<I>})/<I^v  (.)>  (dashed  line); 

(cl  contour  plot  of  the  integrated  energy 
distribution  function  (fraction  of  energy 
greater  than  a  given  intensity  I / 1 . ' , 

IEDF(I/I  .n/'nj;  and  (d)  time  averaged 
distribution  <Ifx.s|>. 


"The  Effects  of  Optical  Smoothing  Techniques 
on  Fi lamentation  in  Laser  Plasmas" 


DISTRIBUTION  LIST 


Dr.  Y.  KatO 

Dr.  K.  Imasaki 

Institue  Laser  Engineering, 

Osaka  University 
Osaka  565,  Japan 

Dr.  John  E.  Murray 
Dr.  John  Lindl 
Dr.  Eric  Storm 

Lawrence  Livermore  National  Laboratory 
P.O.  Box  5508 
Livermore,  CA  99550 

Dr.  J.P.  Vandevender 
Sandian  National  Laboratories 
Albuquerque,  NM  87185 

S.  Kahalas 
R.  Schriever 
U.S.  Dept,  of  Energy 
Office  of  Inertial  Fusion 
Washington,  D.C.  20545 

Prof.  M ,  Rosenbluth 
Institute  for  Fusion  Studies 
University  of  Texas 
Austin,  TX  78712 

Dr.  John  McLeod 
Dr.  Larry  Blair 
MS  548 

Los  Alamos  National  Laboratory 

P.O.  Box  1663 

Los  Alamos,  NM  87545 

Dr.  William  Kruer 
Dr.  Bruce  Langdon 
Dr.  Barbara  Lasinski 
Dr.  Max  Tabak 
Dr.  Alex  Friedman 
Dr.  J.  Mark 
L-477 

Lawrence  Livermore  National  Laboratory 
P.O.  Box  808 
Livermore,  CA  99550 

David  Whelan 
L-45 

Lawrence  Livermore  National  Laboratory 
P.O.  Box  808 
Livermore,  CA  99550 


Warren  Mori 

1-130  Knudsen  Hall 

UCLA 

Los  Angeles,  CA  90024 

Dr.  R.  Sigel 

Abteilung  Laser-Plasmen 

Max-Planck-Insti tut  fur  Quantenoptic 

Postfach  1513 

Boltzmanstrasse  2 

D-8046 

Garching  bei  Munchen  DB 

Dr.  Roger  Haas 
Applied  Physics  Dept. 

University  of  California 

Davis,  CA  95616 

Prof.  E.  Fabre 

Prof.  Guy  Laval 

Ecole  Polytechnique 

F-91128  Palaiseau 

France 

Dr.  Robert  Dautray 

Centre  D'etudes  de  liraeil 

Boite  Postale  No.  27 

94190  Villeneuve,  St.  Georges 

France 

Dr.  M.  Key 
Dr.  F.  O'Neill 
Dr.  D.J.  Nicholas 
Laser  Division 

Rutherford  Appleton  Laboratory 

Chilton,  Didcot,  Oxon,  0X11  OQX,  United  Kingdom 

Dr.  G.L.  Kulcinski 
Fusion  Technology  Institute 
University  of  Wisconsin 
Madison,  WI  53706 

Dr.  S.  Skupsky 
Dr.  R.S.  Craxton 
Dr.  R.L.  McCrory 
Dr.  W.  Seka 

Laboratory  of  Laser  Energetics 
University  of  Rochester 
Rochester,  NY  14627 

Dr.  V.  Happer 
Physics  Dept. 

Jadwin  Hall 
Princeton  University 
P.O.  Box  708 
Princeton,  NJ  08540 


38 


kJil-M.  *Ai 


Dr.  D.  ?3vict 
Dr.  E.  ? rieman 
Science  Application: 
P.0.  Ben  2351 
La  Jolla,  CA  92038 


:n:l.  Corp. 


Dr.  S.  Coggleshall 
Dr.  Roger  D.  Jones 
Dr.  W.C.  Mead 
MS-E531 

Los  Alamos  National  Laboratory 
Los  Alamos,  NM  87545 

L.  Vahala 
Physics  Dept. 

William  &  Mary  College 
Williamsburg,  VA  23185 

R.  Marchand 

Dept,  of  Electro.  Eng. 

University  of  Alberta 

Edmonton,  Alberta 

Canada 

J.P.  Matte 
INRS-Energie 
C.P.  1020 

Varennes,  QC,  JOL  2PO,  Canada 

Monsier  Tudor  Wyatt  Johnston 
INRS-Energie 
1650  Montee  Ste-Julie 
Varennes,  QC,  JOL  2PO,  Canada 

Dr.  Linda  V.  Powers 
Dr.  R.  Berger 
KMS  Fusion  Inc. 

P.O.  Box  1567 

Ann  Arbor,  MI  48103 

Dr.  R.  Davidson 
Plasma  Fusion  Center 
Mass.  Inst,  of  Tech. 

NW16-202 

167  Albany  St. 

Cambridge,  MA  02139 

Code  1220  1  cp. 

Code  2628  22  cps. 


END 

f  ILmeD 

II  ZZ 

D  T  IC~ 


